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CHAPTER I 


GENERAL INORODUCTION 
1.1 IN-IRODUCTION 

The science of physiology deals with the study of 
functional behaviour of living structures. Most of the 
physiological systems involve movoinent of some fluids in 
response to mechanical forces, for example motility of 
the spermatozoa and ova, propulsion of bacteria and 
protozoa, blood flow through circulatory system of raanmals, 
etc. The characteristics of such systems are generally 
studied by using the principles of fluid dynamics and 
making appropriate mathematical models of the system under ' 
consideration, after taking into accoiant the various known 
and unknown variables and parameters related with this 
system. [Sxirton (1965), Gray (1968), Dowben (1969), Fung 
et al. (1972), Rubinow (l975)]. 

In this thesis the following problems have been 
studied. 

(i) Effects of Peristalsis on the motion of micro-organism 
Mathematical models, to study the effects of peristal- 
sis, longitudinal motion of the wall, cilia motion 
etc. on the movement of microorganisms, have beei 
obtained. The results have been applied to study 
the swimming of spermatozoa and ovtam. The effect 
of viscosity variation of the physiological flxiid is 
also studied in this context. 
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(ii) Peristaltic transport of fluid in intestines 1 Effects 
of gradually varying cross-section and multiple ■ 
constriction in a tube have been studied on the 
peristaltic transport of fluid. 

(iii) Blood flow in arteries I Effects of non-uniform 
geometry and multiple stenosis on the blood flow 
through arteries have been investigated by considering 
the non-Newtonian nature of the blood and its 
consistency variation across the lumen. 

1.2 PERIS Ti^S IS 

Peristalsis is the mechanism by which flxiid is 
transported through a distensible tube when contraiction or 
expansion waves propagate along its length. Peristalsis 
appears to be the mechanism for' fluid transport in many 
physiological situations such as transport of urine through 
ureter, food mixing and chyme movement in intestines, transport 
in bile duct, etc. The study of peristaltic transport 
of fluid is based on the well loriown principles of fluid 
mechanics involving interaction of fluid motion in tubes 
with flexible boundaries . In such investigations an approp- 
riate mathematical model of the physiological systa:n is made 
by beeping in view the nature of the biofluid Ci.e. its 
Newtonian or non-Newtonian character, its behaviour as a two 
phase mixture, its viscosity or consistency variation across 
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the lumen# etc*)# the natuxe of the tube (i.e. its elasticity# 
presence of cilia on the inner wall# nonuniformity of its 
cross-section due to constriction^ etc«) and other processes 
(heat and mass transfer, blending and mixing# steady and 
unsteady behaviour# etc. ) involved. 

Several investigators have contributed to the study 
of peristaltic action in both mechanical and physiological 
situations e.g. Burns and Parkes (1967)# Hanin (1968), 

Barton and Raynor (1968)# Rung and Tih (1968)# Shapiro 
et al* (l969)#Yin and Fung (1969)# Chow (l970)# Zien and 
Ostrach (1970)# Li (1970), Lykoudis and Roos (l970)# 

Weinberg et al. (1971)# iylittra (1971)# Jaffrin and Shapiro 
(l97l)# Tong and Vawter (1972)# Brown# et al. (1977)# 

Wilson# et al. (1979), An elaborate review of the earlier 
literature regarding peristalsis has been provided by 
Jaffrin and Shapiro (l97l) • in particular# Burns and 
Parkes (l967) used perturbation techniques to study the 
peristaltic motion through a channel and a tube. Barton 
and Raynor (l968) studied the peristaltic motion in a 
circular tube by using long and short wave length approximat- 
ions. The fluid mechanics of the Ureter has been studied’ 
by Lykoudis and Roos (1970) and Boyarsky and Labay (1972). 

The interaction of Boiseuille flow on the peristaltic 
motion has been studied by M-lttr a and Prasad (1974) . Gupta 



and Seshadri (1976) have investigated the peristaltic 
pumping in a non-tiniform tube. Peristaltic motion of a 
solid fluid mixture has been studied by Kaimal (1977). 

In these studies, the effect of viscosity variation 
of the fluid has not been t^en into account, though 
there have been suggestions that peristaltic mechanian 
may be involved in vasomotion of small blood vessels 
[ Pung and Yih (1968)] ^ in ducts efferentes of the male 
reproductive tract [ Lardner and Shade (1972)] | in 
transport of spermatozoa in cervical canal [ Shukla et al« 
(l978)] and in intestines- all being cases where viscosity 
of the fluid near the wall is different than its viscosity 
in the centre of the duct^ Keeping this in view, Shukla, 
et al. (198 o 4} studied the effects of peripheral layerr 
viscosity on the peristaltic- transport of biofluid and 
compared the results with experimental data obtained for 
flow rates in intes-tine and ductus efferentes of the male 
reproductive tract. 

Though it is known by experiments of Patel et al. (l973'0 
on human feces that it behaves rheologically as non-Newtonian 
power law fluid, little attention has been paid to the study 
of peristaltic transport of such fluids [Raju, et al.(l9T2)]. 
Plcologlou et al. (1973) have noted that peristalsis in 
colon, though rare, is a major propulsive mechanism and 
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presented an analysis to study the movement of hiiman feces 
by considering it as a non-Newtonian power- law fluid. 
Radhahrishnamcharya (1982) has also analysed the peristaltic 
transport of a power law fluid. In view of the power-law 
non-Newtonian behaviour of the chyme and the presence of 
peripheral mucus layer in intestines, Shukla, et al« (1982) 
have studied the peristaltic transport of a power law 
fluid whose consistency varies along the radial direction 
and predicted the effects of peripheral layer consistency 
and the pseudoplastic nature of the biofluid on its flow 
characteristics. 

The effects of peristalsis with particular reference 
to spermatic fluid transport in vas deferens has been 
studied [ Gupta and Seshadri (1976) and Srivastava et al-«(l983)] 
by considering the duct with cross-section varying 
linearly along the length in the stationary frame. 

It is evident, thus, that peristalsis plays an important 
role in the transport of physiological fluids. It is, 
therefore, interesting to study the effects of peristalsis 
on bio-systems like swimming of micro-organisms such as 
spermatozoa and ovum, and on transport of contents of a tube 
with varying cross-section such as intestines with constrict- 
ion due to polyp. Keeping this in view, various models are 
presented in chapters II-V related to these aspects. 
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1,3 EFFECTS OF EERISTALSIS ON THE MOTION OF MICRO-ORGANISMS 

Micro-organi^s like spermatozoa/ nematodes, bacteria 
etc* are known to swim through a fluid by sending lateral 
waves of displacement along their tails or flagella. In 
order for a micro-organism to propel itself with a constant 
velocity, it should move in such a way that the forces 
due to pressure and viscous stresses, exerted by the surround- 
ing fluid on its body, must be balanced. It is also to be 
ixjted that for swimming of these minute bodies, the Reynolds 
number is of order 10“ or less [ Taylor (l95l)] so that 
the propulsion is due predominantly to viscous forces. In 
such low Reynolds number problems the equations governing 
the fluid motion are the linearized Navier-Stokes equations, 
(Stokes equation) fShen et al. (l975)]. 

As the model proposed in this thesis regarding the 
motion of micro-organisms have been particularly applied 
to understand the transport of spermatozoa and ovum in 
hxman genital organs, it would be of some interest to discuss 
the structures of spermatozoon and ovum and related physiology, 
before reviewing the relevent current literature. The proper 
understanding of these phenomena is necessary as it may 
provide an insight into the process of fertility and may 
suggest a biochemical or bicmechanical means to control the 
human population. 
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Structure of Spermatozoon 

The S'permatozoon is a specialised cell which contains 
only the half of the usual 46 chrcanosomes present in other 
cell of human body. Since the nucleus of the female 
gem cell ''ovum' also contain 23 chromosomes, the eventual 
combination of spematozoon and ovum reestablishes’ a normal 
cells of 46 chrcmoscxnes -sdiich is capable of growing by 
division. The total length of spermatozoon is 60 /i-m 
(l Am = 10 ^ meter) approximately equal to the diameter 
of the nucleus of ovum. The structure of a spermatozoon 
is shown in fig, 1.1, Its body can be divided into head' 
and tail. On the basis of slight differences in the 
thickness the tail is fuirther- divided into neck, middle 
piece, main piece and end piece [Warwick and Williams (1973)/ 
Vander et al. (1977), Chatter jee (l980)]. 

The head is oval on surface view, and measures about 
4-5 Am in length and 2«'^-3P,5 Am in dianeter. It is a 
condensed nucleus, containing about 40 o/o DNA bearing the 
sperm's genetic information, and is covered by a cap called 
acrosome. This cap contains several lytic enzymes to enable 
the sperm to pierce the zona pelluicida of the ovum. 

The neck connects head with the middle piece and its 
axial filanents consist of three bundles of fibrils being 
connected with the three basal granules in the anterior 
centriole. 
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The middle piece is the engine room of the spermatozoon 
which controls the whiplike action of the tail» It is 
cylindrical in shape, 5-7 jJ-m long and 1 M-m. wide at base. 

The main piece of the tail is the largest part 
of the spermatazoon. It is 45 Mm long and about 0,5 Mm 
thick at base, gradually tapering towards the end piece. 

It comprises a group of actonysin-like contractile fila:iients, 
arranged in nine-plus-two structures i.e, nine pairs of 
fibrils form a regular array in an outer ring of cross- 
section and two single fibrils are in the center. Contracticm 
of these filaments produces a whip-like movement of the 
tail capable of propelling the sperm. 

The end piece of the tail or Flagellim is the 
terminal portion the axial filament and is about 5 Mm long. 
Its cross-section is identical to that of a cilium. 

MALE REHIODUCTIVE PffifSIOLOGY 

The male reproductive system essentially consists 
of spermatozoa (described above), testes and the male 
genital tract • The testis whose cross section is shown 
in fig.l»2jis primarily composed of many highly coiled 
seminiferous tubules, whose combined length is approximately 
250 meters. These tubules are lined with nimber of cells 
which are continuously regenerating themselves and thus 
paroviding a continuous source of spermatids which are 
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FIG. 1-2 DfAGRAMMATlC CROSS SECTION OF 
A HUMAN TESTIS. 
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ultimately matured as spermatozoa. The other component 
of testis consists of interstitial cells which secrete 
the male sex harmone ^testosterone', necessary for sperm 
production. 

The anatomical organization of the male genital 
duct system is shown in Figure 1,3. From the seminifenous 
tubules, the sperm pass through a network of interconnected 
highly coiled ducts joining to form 'epididymis' vhich leads 
to large thick walled ' yas deferens'. The sperm themselves 
are non-motile at this stage and their transport is achieved 
by means of peristaltic action of contractile cells in the 
duct walls. The epididymis and initial portion of vas 
deferens store spermatozoa, prior to ejaculation. At 
ejaculation the spermatozoa are expelled, by strong contractions 
of smooth muscles lining these duct walls, to be deposited 
at the mouth of cervix in female genital tract. 

STRUCTURE OF OVUM! 

The structure of an ovum is shown in fig. 1.4. It 
consists of a nucleus (germinal vesicle) which is approximately 
equal to 60 ju-m in diameter (the length of a spermatazoon is 
55-60 jJ-m) . This nucleus alongwith the ooplasm centroscme 
is suspended into a sac filled with cytoplasm. This sac is 
surrounded by the zona pellucida and a region "peri vitelline 
space" is *left in between. 
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FtG.1.3 ANATOMICAL ORGANIZATION OF THE MALE 
REPRODUCTIVE TRACT 





FEMALE REIRODUCTIVE PHYSIOLOGY 


The fenale reproductive organs consist of I 

(i) Ovaries ; which are the female counterpart of the 
male testis. It produces ova and fonale sex harmones 
like estrogen and progesterone. 

(ii) Oviduct or fallopian tube I It is a trumpet shaped’ 
duct and is 10 cm. long. The process of fertilization 
takes place in this tube. It is known to have fingerlike 
projections at the ovarian end which sweep oveJ^' the 
ovary during ovulation. 

(iii) Uterus I It is a pear shaped organ. It is on 
average 7, 5" cm long. The length of the uterine 
cavity measures 6 cm. frcan the external os to the 
fundus. Its function is to provide nutrition to the 
fertilized ovum which is now known as blastocyte 
and is grown as embrayo. 

Other sexual organs consist of vagina which is 
approximately 7.5 cm. long and the cervix-which 
connects vagina with the uterus, (see fig. 1.5). 

OVULATION ; 

Each ovary contains a vast number of descrete cell 
clustxires known as primordial follicles. Each follicle is 
composed of one ovum surrounded by a single laij^r of flattened 
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cells. At birth, normal human ovaries contain an estimated 
number of 400,000 such follicles and no new follicle 
is made afterwards. Some 400 of these reach full maturity 
during active sexual life of a woman and all other degenerate. 
The development of a follicle is characterized by an 
increase in size of the ovum and a proliferation of the 
surrounding follicle cells. The ovum becomes separated 
from the follicle cells by a thick mombrane called zona 
pellucida. The ovum is surrounded by a fluid-filled space 
antrvnm which starts increasing enormously and thus pushing 
the ovxm, surrounded by its zona pellucida, towards the 
thin membrane of the follicle . The thin wall of the 
follicle balloons out on the surface of ovary and the 
ovulation occurs when this membrane ruptures. Due to 
the pressure exerted by the antral fluid, the ovum alongwith 
the zona pellucida and tightly adhering follicle cells 
(known as corona radiate) is extruded into the ovarian 
end of the oviduct. The process of ovulation takes 
approximately 2 weeks after the menstrual flow of the 
previous menstrual cyclev 

IHANSPORT OF OVUM AND SP3RMAD0Z0A 1 

At the time of ovulation the ovum is extruded onto 
the stir face of the ovary near the ovarian end of the 
oviduct and is sucked in the oviduct because of the 



17 

sweeping motion of the fingerlike projections at this 
end of oviduct, caused by the contraction and relaxation 
of the smooth muscles surrounding them. The epithelium of 
the oviduct is lined with cilia-hairlike projections 
of the walls, beating in waves towards the interior of the 
duct. This motion of cilia felicitate the motion of the ovum. 
It is speculated that it takes several days for ovum to 
reach the uterus, hence the process of fertilization must 
occur in the oviduct itself. 

The spermatozoa are deposited , during the coitus, 
at the mouth of the cervix and are known to travel at an 
average speed of 1-3 min. /minute. The dimension of the 
fanale genital tract, traversed by the spermatozoa are 
given below 

Length of vagina 7,5 cm. 

Length of the uterine cavity 6 on. 

(from external os to the 
fundus ) 

Length of the uterus 7,5 an. 

Length of the oviduct 10 cm. 

The transport of the sperm to the site of fertilization 
(in the oviduct) is so rapid that the first sperm arrives 
within 45 minutes. This is far too rapid to be accounted 
for by the speirm's own motility. It is suggested that 
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the movement produced by the wave motion along the sperm' s 
tail is probably essential only for the final stages of 
approach and penetration of the ovum. The act of coitus 
provide some impetus because of the ejaculated fluid pressure 
and the pumping action of the penis. The actual mechanism 
of spermatozoa transport is still under speculation. 

The study of swimming of micro-organisms was initiated 
by Taylor (l95l) who modelled it as motion of a fluid 
near a sheet consisting of spermatozoa# down which waves 
of lateral displacement are propagated. Hancock (l953) 
studied propulsion of a thin circular filament through a 
cyclindrical tube# for several wave amplitudes and radii 
of of tube. Further studies were conducted for small wave 
anplitudes and long wavelength approximation by Gray and 
Hancock (1955)# Reynolds (1965), Shack and Lardner (1974), 
Shack et al. (l974), Shen et al. (1975). Attanpts were also 
made to explain the motion of spermatozoa in the fonale 
genital tract by considering the dynamical interaction of 
the wall [Smelser et al. (1973), Shukla et al. (l978)]. Blake 
et al (1983 ) also presented a theoretical model of ovum- 
transport in oviduct incorporating transport mechanisms 
due to ciliary and muscular activity by adding a force 
distribution term in the equation of motion. The motion 
due to ciliary and muscular activity can also be represented 



by a combination of transverse and longitudinal waves 
along the wall. It is also sugges ted ' that peristalsis 
plays an important role in transport of spermatic fluid 
through male genital tract [Guha et al. (1975)# Gupta 
and Seshadri (1976) and Srivastava et al. (1983)]. 

In above mentioned studies the interaction of the 
wall motion on the swimming of micro-organism has not 
been considered. Therefore in this thesis mathannatical 
models are presented to study the motion of self-propelling 
micro-organism when the wall of the duct is undergoing 
peristaltic and longitudinal motion. It may be pointed out 
that the peristaltic motion of the wall could arise due 
to contraction and relaxation of the smooth muscles of 
the uterus and the oviduct while the longitudinal motion 
could be due to cilia on the walls of the oviduct'. The 
effects of peristalsis/longitudinal motility of the wall 
have been investigated by many researchers in other physical/ 
biological situations, f Barton and Raynor (1968)# Fung 
and Yih (l968), Shapiro et al- (l969), Zien and Ostrach 
(1970)^ Shack and Lardner (1972)# Macagno et al. (1975)# 
Shukla et al- (l980ci|jBlake et al. (1983 )] and here we 
investigate these effects to study the motion of micro- 
organism0. As one of the aims of proposed study in this 
thesis is to apply the analysis to fertility control by 
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reducing the speed of micro-organism, the direction 
of the peristaltic waves is taken to be opposite to 
that of its motion vhile studying the movenent of 
spermatozoa and ovtm. 

EFFECT, OF VISCOSITITY OF THE PHYSIOLOGICAL FLUID 

Another important factor, 'affecting the speimi- 
transport and ovum - movement is the viscosity variation 
of the bio-fluid across the cross-section of the duct. 

It is observed that the viscosity of the cervical mucus 
is higher near the walls in the mid - period of the 
menstrual cycle due to concentration of the micelles 
and this aspect should be taken into account for studying 
the motion of the spermatozoa* A detailed account of 
this topic has been given by Scott Blait Cl974), Cliff 
et al« (1950) studied the data related to humans and 
observed the consistency of the mucus at the time of 
ovulation to be minimum . Chemically this change is 
due to the effects of the harmones on the mucus composition. 
The physical aspects of this change have been studied 
for Swimming of spermatozoa in semen and in cervical 
mucus, Lamar et al. Ci 940 > experimentally showed that 
the penetration of the cervical mucus by the spermatozoa 
was maximal for the mucus collected at about the 14 th 
day of the cycle which coincides approximately with the 



time of ovulation in most of the women. This behaviour 
of the cervical mucus was suggested to be utilised as a 
means of fertility control in hxmans [Janes and Marriott (1983)] 

The cervical mucus is a suspension of macromolecules 
in a waterlike liquid (odeblad (1959,1962)]. It is observed 
that in luteal phase the mucus resanbles like a close 
mesh, with macromolecules spreading all over the cervical 
canal with a spacing of 0.03 f.tm between than. However in 
the mid cycle period the molecules group together to form 
micelles which align thanselves along the walls of the 
cervical canal [Odeblad (1968), Elstein (1971), Davajan 
et al. (l97l)]. Because of this alignment of micelles 
along the wall of the duct, a variation in the viscosity 
of the mucus is introduced and its penetration by the 
sperm becomes possible near the axis of the duct. 

To summarize it may be noted that the spermatozoa 
swim by sending lateral waves along its tail while the 
ovum does not have its own motility. Further their motions 
are affected by the contraction and relaxation of the smooth 
muscles of uterus and oviduct. It is also pointed out 
that the female reproductive organs, especially the fallopian 
tube are lined with hair-like projections called ' cilia' , 
whose motion affects the movement of the ovum. This motion 
of cillia is generally represented by a combination of 




transverse and longitudinal waves. Another important 
factor affecting transport of these micro-organisms 
is the viscosity variation of the physiological fluid 
across the duct. In view of these, new models are proposed 
in this thesis, to study, in particTalar, the effect 
of following factors on the movement of spermatozoa in 
cervical canal and ov:3m in fallopian tube. (Chapters II, III, 
IV) : 

(i) Peristaltic motion of the walls 

(ii) Longitudinal motion of the walls 

(iii) Viscosity variation of the physiological fluid. 

1,4 INTERACTION OP NON-UNIPORM GEOMETRY CONSOBICTIONS WITH 
THE PERISTALTIC MOTION 

The second problem studied in this thesis is 
associated with the transport of feces in intestines 
due to peristalsis. It is observed that cross-section of 
the intestinal duct gradually varies along its length. 
Further there may be constrictions due to polyp etc. 
in the limen of the intestine. The dyn^ical interaction 
between the passage of the progressive peristaltic waves 
and an abnormal growth in the Ixmen of the transporting 
duct may be very complex and its study may lead to a better 
understanding of such pathological situations. 
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As stated earlier (in section 1.2), several authors 
have contributed to the study of peristaltic transport 
in physiological as well as mechanical situations' [Burns 
and Parkes (1967) , Fung and Yih (1963), Barton and 
Raynor (1968), Shapiro et al. (l969), Li ^1970), Zien and 
Ostrach (l970), Raju and Devanathan (l972), Tong and 
■Vawter (1973), Manton (1975), Kaimal (1978), Shukla et al.(l978), 
R»K. Aeharya(l982) , Shukla and. Gupta (1982)]. 

It must be observed that in these studies the solutions 
have been obtained in the wave frame, where the flux is 
steady because of the uniform geometry of the transporting 
ducts. However, it is noted that cross section of various 
ducts in human body gradually varies along their length and 
therefore a steady flow through them can not be obtained 
even in the wave frame. Guha et al. (l975) studied the 
effect of peristalsis on the flow through a duct of linearly 
varying cross — section with particular reference to spermatic 
fluid transport in vas deferens. Later Gupta and Seshadri 
(197 6) and Srivastava et al. (l983) also studied this problam 
and obtained the solution in the stationary frame. In these 
studies the value of pressure difference across a wavelength 
has been calculated for a given rate of flow and periodic 
nature of the pressure difference has been shown. 


It may also be pointed out here that the peristaltic 
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transport of fluid in constricted tubes has not been 
studied so far, though such situations may arise due to 
polypt etc, in the lumen of the Intestine, A model is, 
therefore, developed in chapter V, to study the effect 
of non-uniform gecanetry and multiple constriction on 
peristaltic transport of fluids throug^ a circxilar duct, 
with reference to intestines. Following Gupta and 
Seshadri (1976), solutions are obtained in the fixed 
laboratory frame using long wavelength approximation, 

1.5 FLOW OF NON-NEWTONIAN FLUID THROUGH A TUBE WITH MULTIPLE 
STENOSES 

Another important problem dealt in the thesis is 
related to blood flow in arteries with multiple stenoses. 
Stenosis is the term used for the abnormal and unnatural 
growth in the lumen of the arteries carrying blood frcxn 
heart to various parts of the body. It is caused mainly 
due to intravascular atherosclerotic plaques which develop 
at the arterial wall, thus occluding the lumen of the artery 
partly or fully. It has been suggested that hydrodynanict 
factors such as wall stress, pressure distritution and 
seperation phenomenon play an important role in the development 
and parogression of stenosis £Rod]ciwitcz (1977) and Stehbens 
(1982) ], Hence many, theoretical and experimental studies 
of blood flow in stenotic region under different conditions 
have been conducted, [ Young (1968), Caro et al, (1971 ), 



Neron (l974)/ Deshpande etwalr* (1976)# Me Donald (1979), 
Shidcla et.al# (l980a)# Padmanabhan and Devanathan (1981) 
and Talukder et al* (1976)], Most of these analytical 
studies have investigated the flow characteristics of 
blood in arteries, with mild/hon-critical stenosis by 
considering blood as a Newtonian Fluid the si:irvey for 
which can be found in Young (1979), However, it may 
be pointed out that blood being a suspension of cells 
in plasma# has been fotmd to behave as a non-Newtonian 
fluid at low shears rates in tubes of small diameters* 

In view of this some investigators have considered non- 
Newtonian behaviotur of blood £ Shukla et‘ al» (1980b and c) 
Williams and Javadpour (l980), Devanathan and 
Parvathamma (1983)]* Shukla et al* (l980b) studied the 
flow characteristics for mild stenosis by considering 
a non-Newtonian power law fluid model for blood* They 
have also considered the effect of peripheral layer on 
the flow of blood to take into account the clustering of 
cells in the central zone [ Shxdcla et al. (I980b)]* 

In above mentioned studies only the effect of 
single stenosis has been considered and the tube was 
taken to be of uniferm cross-section. However it may 
be noted that many of the blood vessels change their 
cross-section slowly alc^g their lengths [Schneck et al. 
( 1975)3 and may have multiple stenoses at jxanctions, bends 



etc. Talukder et al. (1976) have conducted experiments 
in the case of blood flow in artery with multiple stenoses 
and pointed out that the pressxare drop increases linearly 
with numbers of stenoses. Keeping this in view, in 
chapter VT a model has been proposed to study the effect 
of non-Newtonian power law behaviour of blood and its 
consistency variation near the wall on the flow through 
the artery with varying cross-section having two stenoses, 
one in the uniform region and the other located at the 
junction. 

1.6 SUMMARY 

The thesis is divided into six chapters with 
chapters III and IV being further divided into two 
sub-sections , 

In chapter I, physiological details and brief survey 
of the probl^s and summary of the results are given. 

In chapter II, a math^atical model is presented 
to study the motion of the spermatozoa in cervical 
canal by considering the transverse waves along its tail 
and the transverse and longitudinal motions of the cervical 
wall* The condition of eqxiilibrium of forces on its body 
has been considered. As an attempt to control fertility 
by reducing the speed of sperm, the transverse waves have 
been considered in the direction opposite to the motion 
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of the spermatozoa. It has been shown that by having 
appropriate transverse wave motion and longitudinal velocity 
the sperm may not be able to move towards oviduct even 
if it could continue to have its own propelling velocity. 

In chapter III the motion of micro-organisms of finite 
thickness, in a fluid flowing through a channel and a tube 
is studied under peristaltic and longitudinal motion of the 
walls. The symmetric motion of the transverse waves along 
the surface of the micro-organism is also taken into account. 
The condition of equilibrium of forces on the micro-organism 
is used in making the mathematical model. It has been shown 
that even in the presence of ‘ its own motion, it may not be 
able to move due to presence of transverse and longitudinal 
motion of the wall, when these are applied appropriately* 

This model is used to study the motion of ovum in the fallopian 
tube as a particular case. As the motion of ovim is enhanced' 
due to the presence of cilia motion, it has been shown that 
for fertility control the peristaltic and longitudinal motions 
should be as minimum as possible, i.e. cilia should become 
either immobile or reverse their motion- 

In chapter IV the ssne problem has been studied where 
the effect of the variation of the viscosity due to menstrual 
cycle is investigated* Apart from the above mentioned results 
in chapter III, it has been shewn that the velocity of the 



micro-organiim decreases as the viscosity of the peripheral 
layer increases. 

A mathematical model for peristaltic transport of 
a Newtonian fluid through a txabe with non-uniform cross- 
section and multiple constrictions is proposed in chapter V, 
Ihe shear stress at the constriction is observed to rise 
considerably because of peristalsis and convergence of the 
tube* 

Finally in chapter VI a model of blood flow as power 
law fluid, through an artery having multiple stenoses is 
studied. Effects of peripheral layer and the varying cross- 
section have been investigated. It is found that the shear 
stress and the resistance to flow increase due to increase 
in consistency ef peripheral layer, convergence of the tube 
and the presence of stenoses. 
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CHAPTER II 


EFFECT. OP PERISTALTIC AND LONGITUDINAL WAVE MOTION OP 
THE CHANNEL WALL ON MOVEMENT OF MICRO-ORGANISMS : 
APPLICATION TO SPERMATOZOA IRANSPDRT 

2.1 iniroduction: 

In recent years many investigators have studied the motion 
of self-propelling micro-organisms swimming in a fluid [Taylor 
(1951), Hancock (1953), Gray and Hancock (1955). Reynolds (1965), 
Shack and Lardner (1974), Shack et al. (1974), Smelser et al. 
(1974), Lighthill (1974), (1975), (1976), Rubinow (1976) ]. 

In particular, Taylor (1951 ) studied the propulsion of spermat- 
ozoa by modelling it as a two-dimensional, infinite extensible 
sheet of zero thickness with a sinusoidal wave travelling down 
its length/ and Hancock (1953") considered the self propulsion 
of a thin circular filanent and calculated its projxilsive 
velocity. Later Reynolds (1965) studied the swimming of micro- 
organism in a fluid for small Reynolds number and small wave 
amplitudes. Shack and Lardner (1974) also studied the motion 
of a micro-organism represoited by a thin sheet, in a channel 
under long wavelength approximation and have shown that its- 
speed increases as the amplitude of the propelling wave increases. 

In above mentioned studies the interaction of the wall 
motion on the swimming micro-organism has not been considered. 
Therefore, in this chapter, a mathematical model is presented 
bo study the motion of self-propelling micro -organi an when the 
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wall of the channel is undergoing peristaltic and longitudinal 
wave motion. It may be pointed out here that the peristaltic 
motion of the wall could arise due to contraction and relaxation 
of the smooth muscles of the uterus and the oviduct while 
the longitudinal motion could be due to the presence of 'cilia', 
the hair-like projections on the walls of the oviduct. These- 
effects have been investigated by many researchers in different 
situations [Barton and Rayrxjr (1968), Fung and Yih (1968), 

Shapiro et al. (1969), Zien and Ostrach (1970), Shukla et al» 
(1980), Shack and Lardner (1972), Macagno et al. (1975)]. 

One of the aims of the proposed study is to apply the 
analysis to fertility control by reducing the speed of spermat- 
ozoa in female genital tract. It is therefore suggested that if 

^ « 
the peristaltic waves are induced in the direction opposite to 

motion of the spermatozoa by some biochemical or other means, 
it may be possible to reduce the speed of spermatozoa to such 
an extent that it may not be able to reach the point of fertiliz- 
ation and thus avoiding fertilization. 

2.2 MATHEMATICAL MODEL 1 

Consider the swimming of a thin propelling sheet in 
a Newtonian incompressible fluid flowing through a two-dimensional 
channel having flexible boundaries. It is assumed that the 
sheet, while swimming, sends down lateral waves of finite 
amplitude along its length. Furtlr^r, peristaltic waves of finite- 




Fig.2*1{a) Propagation of an elastic sheet 
swimming through a two- 
dimensional channel with peristaltic 
and longitudinal motion of walls 
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amplitude are imposed along the flexible walls of the channel 
in the direction opposite to the motion of the sheet as shDvjn. 
in Fig. 2.1. The sheet is considered to be swimming with a 
propulsive velocity in the negative axial direction. It 
is assumed that the waves travelling along the channel walls 
and along the sheet are in synchronization under steady state 
and thus have the* sane wave speed C (along positive axial 
direction) and the wavelength In a fixed frame of reference 

(X',Y',t^) the shapes of the walls and the sheet at an 
instant t* are given by, 

H'(X",t') = at b,sin~(X'-Ct'+V' t' ) (2.l) 

1 1 K p 

= -a t b^sin~(X'-Ct'+V' tO (2.2) 

^ ZAP 

H' (X' , t' ) = b^sin|^(X' -Ct' tV* t' ) (2.3) 

o O A p 

Where b-]_, ^^2 amplitudes of the peristaltic 

waves travelling along the upper wall and the lower wall 
respectively, 2a is the width of the channel and b^ is the 
amplitude of the lateral wave through the sheet. 

Since the Reynolds number involved in the swimming of 

—3 . 

a micro-organism is of order 10 , the governing equations 

for the fluid flow can be written, after neglecting the inertia 
terms, as follows 

. ir + iu' ? "y; ] = o 

9x' dx' dY' 


(2.4) 
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3y' 3x' ^ 


^ 0 


(2.5") 

( 2 . 6 ) 


3X' 3Y' 

where are the vel9cit;y components along the X* and Y* 

directions respectively, P' is the pressure and 4 is the 
viscosity of the fluid. 

Taking into account the longitudinal motility of ther- 

# 

^^ralls ,the boundary conditions are I 


U' (x' , t' ) = 

Gi(X' 


an' 

V' (X' , t' ) = 

1 


at' 

U' (X', t' ) = 

G^(X‘ 


aHi 

V' (x' , t' ) = 

2 


3t' 

u"(x',t') = 

-V' 

■ P. 


dH' 

V' (X" , t' ) = 

_p 


at' 

where G£ and G' J 


P 


at Y' = H' (XSt') (2.7) 


[. at Y' = H' (X' , t' ) 


( 2 . 8 ) 


at Y' = (X' , t' ) 


(2.9) 


of the walls. 

Further the sheet is self-propelling and the forces 
exerted, by the fluid, on it must balance for its motion with 
a constant velocity, i.e. 

/ (x"*'+T )tlS =0 (2.10) 



42 


mm. 

where t and r are the resultant of forces acting on the 
upper and lower surfaces of the sheet respectively. 

In a frame (x'/y'/t') moving with velocity C-V' in the 
positive axial direction, the sheet and the walls of the 
channel appear stationary and the flow in this moving frame 
steady. 

Using the following, transformations 


x' = X'-Ct'+V't' ; y' = Y' ( 2 . 11 ) 

P 

and introducing the following dimensionless quantities 

X = x'/A, y = y'/a, t = Ct'/X, p = Ef/ ^^2" 

a 

(u,g3_,g2#Vp) = (u',g£,g^,v^)/c, v = ^ 


(where symbols denoted by small letters in the moving frame 
represent the same quantities as denoted by capital letters 
in the fixed frame) we can write dimensionless form of eqis. 
(2. 4 ’) -(2. 6 ) in the moving frame aa follows 


. |E + p 2 ^ 0 

3x^ 3y^ 


^ + p2(p2 tv AY) = 0 


ay 


ax 


dy 


^ + -AY 
3x By 


0 


( 2 . 12 ) 

(2.13) 

(2^4) 


where. 


P = a/ A. 


(2.15) 
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The boundary conditions for u are written in 
dimensionless form as 


where 


u 

u 

u 


Vp + g^(x) - 1 
Vp + g^Cx) - 1 

-1 


; y = h^(x) 
# y = h 2 (x) 
,* y = h^(x) 


(2.16) 


h^(x) = sin 211 x 

hj^(x) = 1 + sin 2iEx (2.17) 

h2(x) .= -1 + ^2 2 It X 

For long wavelength approximation (P<< l) the eqns. (2il2) 
and ( 2 . 13 ) (after neglecting the terms involving h and hdgher 
powers) are reduced to 


3P.. 

d X 



6Y 


- P = 
dy 


o 


(2.18) 

(2.19) 


Further using the stress-strain relationship to obtain 
the stresses in terms of u and v and applying long wavelength 
approximation /the force equilibrium condition (2,10) can be 
written as 

1 

/ [p] dx = 0 ( 2 . 20 ) 

o 

^ ^!.fy ]y=h (x) 8x ° ‘2.21) 

o 

where [f]* indicates the difference ift the <^antity f 
above and below the sheet . 




The equation (2.18) and (2.21) are solved in the following 
section, using conditions (2.16), to obtain the propelling 
velocity of the sheet and the total fluid flux in the charnel. 

2.3 MATHEMATICAL ANALYSIS 


The differential equation (2. i8) is solved with boiindary 
conditions (2.17) to obtain the velocities u"*” and u“ in the 
regions ^ ~ ^ < V respectively, which 

are given as follows 


u 





(h, -y) 


-1' -S' + u 


(2.22) 

• V^o ° 

(hi-h^) 


(y-h ) 

(h2-y) 

(2.23) 

h2-:h - ""o 

(h'-h^)' 

2 o 


The dimensionless flow flux q(= q^/ac) in the moving frame 
can be obtained as 

^1 

q = / u dy 


q” = / u" dy 

^2 

which on using eqns. (2.22) and (2-23') give 


q 


*■ = - )(h^-h.,)^ + i(h^-h2) (u^+u^) 


12 'dx 


" 2 - 0 - 
+ 

ftp” 

3x ' 


(2.24) 

(2.25) 


To obtain the pressure gradients 


the eqns. (2.24) 



and ( 2 . 25 ) can be rewritten as 


= -12 - + 6 

(hn-h (h,-h 

X O 1 o 


F" = -12 

dx 


+ 6 


V“2 


(h -h,)^ (h -h,)^ 

o 2 o 2 


(2.26) 


(2-27) 


Integrating equation of continuity (2.14) it can be seen 
that the fv?o fluxes are constants. Further frcxn eqn. (2.19) 
it is evident that the pressure rise over a wavelength Ap is- 
same for the two regions. 

Integrating eqns. (2.26) and (2.27) over a wavelength, 
we get the following two relations among the three unknown 
quantities V , q and q 

6V12 
-12-3"l21 6V22 

Using the expressions (2.22), (2.23-) and (2*26), ( 2 * 2 T) 

for u— and 9 ^ / the force equilibrium condition (2*2l) pro- 

amm 

vide us the third relation, necessary to determine V^,q and q , 


+ 61^3 -I 2 I 32 = Ap 

(2.28) 


( 2 . 29 ) 


as 


I2q'*’l2^ - I2q“l22 + ^ ^^ 3 4 +^ 35 ^- 1212 g = 2Ap (2.30) 

where . 


'll 


^12 


^13 


= / 

0 

_ax 

I 

II 

H 

CM 

H 

1 

,._dx 

* 


= / 

0 

1 

(h,-h^)^ 

1 0 

g., (x) 

i 

^22 ~ 

= / 

— i . . dx 

m ' 

I.... = 

0 

(h,-ho)^ 

f 

23 


0 

1 
; 

0 

1 


dx 

(h -h.,)' 
o 2 

_dx 

(h_ -h, ) ‘ 
o z 

g^ix) 

(h -h^) ' 
o 2 


(2.31) 


dx 



4 b 


and 


31 


1 

= / 
o 


h, +h -1 

-tv» 


dx 


^3 2 

o 


1 

2 O 


(h -h_)3 

O 2 


dx 


3-2h, -4h 
1 o 


'•33 ~ " {2 

o (hT-h } 

1 o 


- -h 


3+2h«-f4h 

__ 2 p 

(h -h 
o 2 


] dx 


1 3-2h,-4h 

1 3+2h„+4h 

'35 " it (f, 'Tja ] 32 

O 2 


1 1-h, ~h 

3® o (h,-h 

1 o 


1+h^+h 

+ -■?-■ °-jax 

(h-h.r 
o 2 


(2,32) 


Solving equations (2. 2.8) -(2 .30) we get 






3 2 

9'i 

I, 


( 6 l 23 ~ 12 I^,)] 


22 


[2I33 + "■ ^22^ 


(2. 33 ) 


and the total flux q = q + q" as 


r 1 L ^ ^ ^ ^22^^ . ^13 ^^12 . ^23 ^^22 v 

^11 -^21 -^11 -^21 ^ ^11 ^21 


(2.34) 


The flux in the stationary frame Q is related with the 



flux q in the moving frame by the following relations 


Q = q + (l-Vp) 

Averaging equation (2.35) over a time period we. 


(2.35) 


therefore, get the time averaged flow Q in terms of q and V as 

ir 


Q = q + 2 (1-Vp) 


(2.36) 


Sqns. (2. 31) -(2. 3 6) provide expressions for ^ » q and 
Q for the given forms of h^ (x) , h^ (x) , h 2 (x) (eqn.2,17) and 
general forms of functions g^(x),g 2 (x) representing the longit- 
udinal motility of the channel. In the following analysis we 
consider the following forms of longitudinal motion of the wall* 


g^ (x) — ^11 ^12 ® 2 X 

g 2 (x) = + C 22 sin 2 'll x 


(2.37) 


For these given forms of h^,h^,h 2 and *32.^92 
writing 


= e,-e , a = e -e 
1 1 o 2 o 2 


(2.38) 

The expressions (2.31) and (2.32) can be rewritten as 


I ^ ; ... ,„„dx , 

o (l+a^sin2 t x)^ 


I = / -- - ' (i=l,2) 

o (l+cc^sin2 n: x) 


i3 xX x2 x2 x4 


o (l+a^sin2 m x)^ 


(2.39) 



and 


31 


3 2 


33 


^34 


35 


36 


■37 


3 8 


(e +E, ) f - -,ax 

o (l+a^sin 2 it x)'^ 

'(e„+e,) / ax 

o (l+a2sin 2 tt x) 

I 12-2 ^£i+2e^)lj_^+l22+2 

^^12“^ ^14“2 ^£3_+2e^) '^2l'^^12 

It ^2 2"^^ ^ ^ ^24^ ^2l'*’f“ ^2 4'*'^ ^ ^ 2*^^ ^ ^3 8 ^*^22 

-(£l+eo)l^4 + ^^2 ‘^^o^^24 

/ -_sini.2 7.„x „ 

0 (1+a^ sin 2 u x) 

1 . 2 „ 

j. . -2 n_x^_. 

o ^i+^2 2ltx)‘^ 


(2.40) 


These integrals can be evaluated by contour integration and 
are listed in appendix 2 (a) . 


In 

by using 


the following, scrae particular cases are obtainetd 

the following relations in the integrals, for integer m,n 


1 

o (l+asin 2 "n: x)^ 


1 

/ 

o 


dx 

(l— asin 2 t x)^ 


and 


1 

/ 

o 


m 

sin_ 2 TE x__ dx_ 
(l+asin 2 ic x) ^ 


(- 1 )™ 


1 

/ 

o 


sin^ ^TCxdx 
(1-asin 2 itx)’^ 


( 2 . 41 ) 


When there is no peristaltic motion Of the walls, we can write 



and using relations (2.4l) we obtain the following relations in 



the integrals given by (2.39) and (2.40) 


^11 “ ^21 


^12 “ ^22 


^13 ~ ^11^12'^^12^14 


T = r T —r* T 

23 ^21 12 22 14 


"l4 = '^24 


^31 “^32 


^33 “ ^^12 " ®^o^l4 


"34 

- Cj_3_(I^2”4£qIj_4) 

^12^"l4"'^^o"37 

"35 

H 

H 

0 

u) 

1 

CM 

H 

H 

H 

CM 

0 

It 

“ S2^"i4“^S"37 

"3 6 

= -2e^I, . 

0 14 


I> 

CO 

H 

~ "38 



(2.42) 


We therefore get the expressions for V^, q and S from equations 
(2.33), (2. 34)^ (2.36; as, 

1 


V- = j;f(M Ap+M) 

o 


(2.43 ) 


, 6M, M , 

= gAp(- M“^I^Y^'*’2I^ ■I^^12^^1l‘^^2l'^^14^^12‘’S2^ 

11 o 11 o 11 11 


“^^^12 ^ 2 " I ^"'^ C ^^ 12 ^‘^ 11 '^'^ 21 ^'^^ 14 ^'^ 12 ”^ 22 ^“‘^^ T 2 ^ 


12 - 


(2.44) 


n _ 1 A f 1 4.(^ i2 ^12 „v 

Q - 6 Ap(- 2 - +6^ ^ - -12^ )+ ^ ( j ' -2) + 2 
11 o 11 o o 11 


■*■ ^1 2 ^^ ll '‘'^2 1 ^ '‘■^14 ^‘" 12“‘^2 2 ^ ^12 ^ 


(2. 45') 
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where 

I 

M =2 + 2 - 

_ 11 


-24^14 i:^12«^ll-^=21^«14‘=12-S2’] 


( 2 . 46 ) 


Further if longitudinal motion is also absent i.e- 


g^(x) = 0 , ^2^^^ “ ® 


we get a case when the sheet is swimming through a channel. 
Using the values of the integrals frcsn appendix 2 (a), we get 
eqns. (2 .43 )- (2 .45 ) reduced to 


, = i(l-e2)V2 Ap . 



2 

l+2e^ 

2 (l-e^) 
o 

2 

l+2e 


(2.47) 


o 

Q = - “(2+e^) (l-£^)^/^ Ap 

These expressions are sane as obtained by Shack and 
Lardner (1974) for q = 2Q and Q = 2Qp. Thus the proposed model 
includes the case studied by Shack and Lardner (1974) as a 
particular case when there is no peristaltic and longitudinal 
motion of the walls of the channel. 


2.4 RESULTS AND DISCUSSIONS 

The values of V and Q have been conputed using eqns. (2.33) 
(2.36) for the following set of parameters 



51 




= 0 , = o, = (0,0,2,-0.2),C^(.=C^^)=.{0,0.2^-.2} 

= (0,0.35), ^2 =C0,0.25>-0,25),Ap =(0, .05, -. 05 ) 

and are plotted in figs 2.2 - 2»13 as functioiBof 

Effects of various parameters on is graphically 
shown in fig. 2.2 - 2.7. In fig. 2.2 the results corresponding, 
to the particular case, when there is no peristaltic motion 
on the walls (e-j^ =0, = o), are shown. As noted in the 

expressions (2.43)~(2.46) the velocity amplitudes C^{=C^ 2 ^ 
and '^2^~'^22^ appear always as their difference C 2 —C^=d. The 
effect of this parameter d is therefore, shown for Ap = O, . 05 . 

It is observed that increases as magnitude of increases. 

For fixed 0, decreases as d increases and for < 0, a 
reverse trend is noticed. It is also noted that, in this 
particular case, there is no effect of the longitudinal velocit- 
ies ^^,^2 If =0 which is also evident frcxn the expressions 
( 2 . 43 ) and ( 2 . 46 ). The effect of Ap, for fixed values of 
and d is to increase V^. 

Effect of and C 2 for non-zero values of e^, £2 is. 
shown in fig. 2.3 for Ap = 0. The effect of and is similar 
as observed in fig: 2.2. Comparing fig. 2.2 and 2.3, it can 
also be noticed that because of the peristalsis on the walls, 
the propagation velocity V is reduced and, in fact, is in the 
opposite direction. The maximum negative value of is obtained 
for the combination < O, O , 
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d = 0'4 



Effect of Qnd C2 on Vp 
ef Sc £2 peristalsis) 


L - 
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feet of Cl and C2. on Vp for fixe 
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/lEr-o 

I £2=0 


„ 4 El = 0-35 
//0-3 i £2=0-25 


AP =-0-05 
ap = 0 
6P = 0-05 



Similar effects of ^2.* *^2 seen fron 

fig» 2.4 - 2 - 6 , where effects of e-^ and £2 shown for fixed' 

(Cj^,C2) = [K 2 /. 2 ), (0,0)^ (, 2 ,-, 2 ) ] and Ap = 0 , Comparing 

these three figures it can be observed that for = *^2 “ 
sign of £2 does not matter. As for the maximum negative value 
of Vp for fixed £q/ S j. concerned, it is obtained for the 
first case CC2^ ^ 0 , C2 ^ 0 ^ with £2 ^ ^ Cfig. 2 -. 4 *) * 

Increase in Ap, the pressure rise causes increase in V . 

ir 

A negative value of Ap, however, reduces the value of V to 

Jr 

negative values (fig. 2 . 7 ). 

Effect of various parameters on time-average flux 5 
is shown in fig. 2 • 8 - 2 . 13 . It can be seen that in general the 
results are opposite to those obtained for V . The effect of 

ir 

peristaltic waves on the wall is to increase D but increase in 
the magnitude of the amplitude of the causes a decrease in Q, 
(fig. 2.8,2.10,2,11). It also increases as Ap becomes negative 
(fig. 2.12,2.13). For a particular case of no peristalsis 
(e^ - ^2 ~ effect of d is oppositie to that of i.e. 

Q is minimum for d> O, (fig. 2 . 8 ). However in general case 
maximum value of S is obtained when and C2 are both negative 

(fig . 2.9). 

2,5 MOTION OF A THH'I SHEET IN A CHANNEL UNDER PERISTALTIC MOTION 
OF THE WALLS 

This case can be obtained from the present model when 
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Fig.243 Effect of Q ( Q =*2 ,C2=-2) 



there are no lateral wave on the sheet and no longitudinal wave 
motion along the channel walls. It should be noted that this 
case is physically different fran the usual peristaltic flow 
in a channel [Shapiro et al^ (1969)] as there exists a sheet 
in the middle of the channel which divide it into two separate 
zones , 


Mathematically, this case is obtained by putting 
g^^Cx) = 0, g 2 Cx) = O and = O, in the equations (2 ,33 )~(2.36) 
and then the expressions for V , q and 0 become (for the symmetric 
case when = -^2 = 






■11 


11 


M,I, 


T; 1 a ^ 1 a 112 

Q = T ApC--f -i-6 


where 


'11 
M = 41 


o 11 

I. 

_ j 

I, 


I, , MI M 

2— ^)+ C-2-' ■ +2H — ^ —2.—^) 

M 1^, -^^M ^ ^ ^I,- ^ M I^T ^ 

o 11 o 11 o 11 o 


(2.4S) 


12 ” ®^^14 


Ml = 2 + 2 I 


3_1 

11 


(2.49) 


«2 = ^14 + 4'^12 

The effects of Ap, and £2 ^ ^ seen 

frcsn figs. 2.2, 2. 5,2. 8 and 2.12 at = 0. It can be noted 

frcxn fig. 2.2 and 2.3' that the velocity of the sheet increases 

as Ap increases. However for £ 2 ^ 0, £2 ^ 0 it is seen that 

Vp is negative, i.e. the sheet flow in the direction of 
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peristaltic waves and this effect is enhanced as and £2 
increase. The effects of and £2 on Q can be studied 

from figs, 2.8 and 2.12. It is observed that Q increases 
as the magnitude of and £2 increase and as Ap decreases. 

2.6 APPLICATION TO SPERMATOZOA TRANSPORT I AN APEROACH TO 
FERTILITY CONTROL 

The spermatozoa are male reproductive cells. They are 
specialized cells containing only 23 out of full 46 chrcmoscsnes 
of a normal cell in human body. During coitus, the spermatozoa 
are deposited, at the time of ejaculation at the mouth of the 
cervix in fanal genital tract. It has been found that tliey take 
on average 45 minute to reach the fertilization point in the 
fallopian tube and thus travel with an average speed of 1-3 m.m./ 
minute. It is necessary, from the point of fertility control, 
to reduce this speed of spermatozoa in the female genital tract 
by some biochCTiical and other means. 

It is known that spermatozoa, when in large number and 
close to each other, travel in unison. It is therefore possible 
to approximate their motion by a sheet flowing through fluid 
[Taylor (1951)], The proposed model and the results discussed 
in section 2.4, can therefore be applied to study the effects 
of peristaltic and longitudinal motion of the walls of the 
fan ale genital tract on the spermatozoa transport. 

Noting the effects of various parameters on V^, as 
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discussed in section 2*4, it is evident that V can be made 

P 

negative by application of a negative Ap (i.e* pressure drop) 
and this effect can further be enhanced by inducing tte 
peristaltic waves travelling along the walls in the direction 
away from the point of fertilization. Further the motion of 
the cilia can also enhance this effect by a suitable choice 
of longitudinal wave. It is therefore suggested that if the 
generation of pressxxre drop and the peristaltic motion is 
possible by some mechanism, such as sudden withdrawl of male 
organ from the vagina immediately after ejaculation, it may be 
possible to reduce the motion of the spermatozoa and hence 
controling fertility* Thus the ancient belief behind the 
practise of withdrawl of the male organ immediately after' 
ejaculation, to control fertility, may have a scientific basis 
as predicted by this model. 

2.7 CONCLUSION : 

A mathematical model to study the effect of peristaltic 
and longitudinal motion of the walls, on the propalsion of a. 
sheet has been presented and the results have been applied 
on the swimming of spermatozoa through fenale genital tract. 

It has been shown that the speed of spermatozoa can be reduced 
considerably, to avoid fertilization, by generating a pressure- 
drop and inducing peristaltic waves on the wall of the genital 
tract through some biomechanical and other means such as 
sudden withdrawl of male organ fronitie vagina, immediately 
after ejaculation. 
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APPENDIX 2(a) 
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CHAPTER III 


MOTION OF MICRO-ORGANISMS UNDER PERISTALSIS :APELICATICN 
TO OVUM IBANSPORT 

3.1 iniroduction: 

In the past few decades the fluid-dynsnical aspects 
of peristaltic pumping have been extensively studied because 
of its importance in transporting fluids in many biological 
systems [ Biirns and Parkes (1967), Bart(Xi and Raynor (1968), 
Fung and Yih (1968), Shapiro et al. (1969), ChDw (1970), 

Jaffrin and Shapiro (l97l), Weinberg et al. (1971)]. A review 
of much of the literature regarding this has been presented 
by Jaffrin and Shapiro (l97l) . Most of these studies pertain 
to the motion of urine in ureter and chyme movement in intest- 
ines and very little attention has been paid to study the 
effect of peristalis on biosystens such as the transport of 
micro-organisms like spermatozoa and ovum. Further it is 
also known that cilia attached to the surface of transporting 
ducts, such as oviduct, play a significant role in transport 
of fluids in many situations [Barton and Raynor (1967), Lardner 
and Shack (1972)]. It would, therefore, be interesting to 
study these effects on the transport of micro-organism such 
as ovum. 

The Study of self-propelling micro-organisms was 
initiated by Taylor (1951 ) who assumed the organism to be 
an infinite sheet of zero thictoess. Later Reynolds (1965) 
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extended this analysis to include the effects of inertia* 
Shack and Lardner (1974) also studied this problon under long 
wavelength approximation. Att^pts were also made to explain 
the motion of spermatozoa in the female genital tract by 
considering the dynamical interaction of the wall £Smelser 
et al . ( 1973 )/ Shukla et al. (1978)]. 

Recently Blake et al* (1983) have presented a 
theoretical model of ovum-transport in oviduct incorporating 
transport mechanisms due to ciliary and muscular activity 
by adding a force distribution term in the equation of motion. 
The motion due to ciliary and muscular activity can also be 
represented as a ccmbination of transverse and longitudinal 
waves along the wall and effects of this interaction should 
be considered on ov\m transport. 

Keeping these in view, in this chapter models are 
proposed to study the effects of peristaltic and longitudinal 
wave motion along the wall of the duct (rectangular and.. 
circular) on the locomotion of micro-organisms of finite 
thickness which may also generate symmetrical peristaltic 
waves along their bodies. The governing equations of this 
syst®n have been solved under long wavel^igth approximation 
and the effects of various parameters on the velocity of 
micro-organism and the flux, are studied. This analysis 
has also been applied po study the transport of ovum in 
oviduct. As one of the alms of the proposed study is to 
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suggest a biomechanical or biochemical means for fertility 
control, the interaction of peristaltic waves and longitudinal 
motion and pressure difference on the movement of ovum in 
the oviduct have been studied. It has been shown that if 
there are no transverse and longitudinal waves along the 
wall of the oviduct, the ovum may not be able to move forward 
for fertilization. 

3-2 TRANS EDRT THROUGH A CHANNEL (CASE l) 

3.2.1 MATHEMATICAL MODEL 

Consider the propulsion of a flexible slab (micro-organia:n) 

of infinite width and finite thiclcness 2$' in the middle of' 

a channel. The space between the boundaries of the slab and 

the channel is filled with a Newtonian fluid and the walls 

of the channel are undergoing longitudinal wave motion. It 

is assumed that the slab is being propelled at a constant 

speed V' in the negative X' -direction. Let travelling 
P 

sinusoidal waves with speed c-V' and wavelength X but different 

ir 

amplitudes b' and b^ travel along the surfaces of the channel 
and the micro-organism, respectively- In a stationary frame 
(x*,Y',t') the shapes of the wall of the channel (Y' = + H' ) 
and the micro-organism (Y' = + H^) at any instant t' are 
therefore given by (see Fig. 3.l) 

H'(X',t') = a + b' sin ^ (x^ - ct' + \^t' ) 


(3.1) 
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IG. 3.1 


y 



\ 

t 



(b) 


SWIMMING OF A MICRO ORGANISM THROUGH A CHANNEL 
(a) WITHOUT PERISTALSIS (b) WITH PERISTALSIS 





(3.2) 


H'(x%t') = 5' + b' sin ^(X' . Ct' + V t') 

i A P 

Since the Reynolds number, as in Chapter II, is of 
order 10 , the equations governing the flow of the fluid 

in the gap between the boundaries of the microorganism 
and channel (after neglecting the inertia terms) are 


i? 6ir , 

ax' ax' ay' 

(3 .3 ) 

^ M i = 0 

(3 .4') 

dY' 6x' ^ 3y' 


= 0 

ax' dY' 

(3.5-) 


where U',V'' are the velocity ccxnponents of the fluid, 
is viscosity and P* is the pressure. Talcing into account; 
the periodic nature of the longitudinal motion of the channel 
walls, the boundary conditions are 

U' (x',tO = G' ix' - ct' + v't) = c, sin ^ (X'- ct' + V' tO , 

(3,6) 

V' (X' , t' ) = + ‘ ; y' = + H' (X' , t' ) 

and U' (X' , t-^ ) = - , 

V'■(X^t') = + ‘ ^ ; Y' = + H'(X',t') 

“■ dt ~ 

Under the present symmetric situation, the force 
equilibrium condition on the surface of the micro-organism 
can be written as 

/ t' 

S 



dS = 6' AP' 


(3.8) 



where t * is the resultant of the forces acting on the 
surface of the micro-organism AP' is the pressure rise 
over a wavelength and S is the surface of the micro-organism. 

In a frame moving with velocity c-V' in the positive 
axial direction, the boundaries of the channel and the 
micro-organism appear stationary. Using the following 
relations to transform various quantities fran the stationary 
frame Cx',Y',t') to the corresponding quantities in the 
moving frame (x, y, t) in the dimensicnless form, 

X = (X' - ct' + V^t')/X, y = Y'/a, t = |t' 

u = (U'- C- v;)/c, V = VV ¥ , P = ^3.9) 

P ^ 

(g, c^, Vp)=(G^,c',V^)/c, (h,h^, 6 ,e,ej^) = (h' ,h', 6 '/b',b' )/a 


eqns 


. (3,3) - (3 ♦ 5 ) are transformed as. 


. .aP. + in ^ Jn ^ 0 
dy? 3y^ 

. 8B-+ a2(p2 d^Y) ^ 


dx 


0 


6 x 


ay‘ 


(3.10) 
(3 -.11) 


+ . 3 Y- = 0 (3.12) 

6 x 3 y 

where P = a/X, 

For large wavelength (P << 1 ) we can neglect the terms 
involving P^ and higher orders and thus the eqns, (3.IO) and 
(3,11) are simplified to. 


(3-i3) 
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- = 0 
3y 


(3.14) 


The boundary conditions for u are obtained from (3.6) 
and (3 . 7 ) as 


^ = Vp+g(x)-l ;y = + h(x) 


( 3 . 1 *^) 


u = -1 


^ y = + hj^(x) 


where 


h(x) = 1 + e sin 2 11 X 
h^(x)= 6 + sin 2 TC X 


(3.16) 


g(x) = sin 2 15 X 

Using the stress - strain relationship and neglecting 

terms containing and higher order, we obtain from (3.8) 

1 

(3.17) 


/ p dx = 0 

o 


dh,(x) 


and ^ = h^U) + dx‘ 

Using the symmetry we solve these equations fon- 
y ^ O region only. 

3.2.2 ANALYSIS. I 

The equation (3.13) can be solved, tising boundary 
conditions (3.15"), as 

U = |( -(w-hj^)y + hh^] + - 


(3 . 18 ) 


1 


( 3 . 19 ) 
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The volxjmetric rate of flow is 
h 

q = / u dy = - •) (h-h^)^ + |(g + Vp-2) (h-h^) 

(3.20) 

The expression (3. 20) can be rewritten as 


— —12 

ax ~ 


, q . . ^ 

(h-h^)^ (h-h^)^ 


( 3 . 21 ) 


Integrating equation of continuity (3.12), it can be 
-aeen that the flux q is constant in the wave-frame. Integrating 
equation (^.2l) over a wavelength we get 

Ap = -12I^3_q + 61^2'^P + *■ ^21^2 (3.22) 

dx ^ \ dx ^ - V 


^ dx ^ 

where 1.1=/ , I,. = j* , I,. = ^ ^ ^ dx 

o (h-hi)3 ^ (h-h^)2 13 o (h-h^)2 


(3.23) 


Using expressions for u and from (3.19) and (3.21) the 
force equilibrium condition (3.18) can be written as 


^^21^ " ^^22^p “ ^^23 ''' ^^24 


h+h. 


1 ht2h. 


where I 


= / 


,dx , 


21 “ ' ^ ' *22 
o (h-h^) 


= I 


o (h-hj^)' 
1 h+h. 


(3-24) 


(3 .25) 


‘23 


1 (h+2h,)g(x) 

= r — ' o dx - I^ - = / - - 

„ (h-h,)^ o Ch-h^)^ 


The equations (3.22) and (3.24) are linear equations 

in V and q and can be solved as 
P 

Ap. + 


'*11 
oo 


(3.26) 



q = 

00 



(3 

.27) 

^00 

^^12*^21 ** 

^U2‘^U 




^11 

= "21 





^12 

= -6l2;j_.I^3 

+ + 


^24 


^21 

3 ^22 



(3 

.28) 

^22 

= 2112^^23 - 

3 I 24 ) + 21^3 

^^^12 " ^13^ 




The fluxes in stationary and wave frames are related as 

(3.29). 

The time - averaged flux in stationary frame Q / therefore 


Q = q _ 


is related with q. and as 
Q = q - (V -l) (l-<> 5 ) 


(3 .30) 


The integrals (3-23) and (3-25) can easily be evaluated' 
in terms of the integrals given in appendix 2 (a) (Chapter II) 
and using eqns. (3*26) — ( 3 . 30 ) and Q can be obtained. 

3.2.3 RESULTS AND DISCUSSIONS I 

The expressions (3*26) - ( 3 . 30 ) are ccmputed for 

various values of c^/e,Ap,e ^,6 and the effects of various 

parameters on V and Q are shown in fig. 3*3 — 3.11* Effects 
P 

of c, ,e Ap on V as function of are shown in figs. 3.3 - 3.7. 
It can be noticed from these figures that 


V 






U A. 



FIG. 3-4 


effect of and t ON 


Vp (Cv=0/S = 0.2) (CHANNEL) 
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£=-.2 



-15 


FiG. 3-9 EFFECT OF Ap AND € ON Q (CHANNEL) C C| =0.6 =0*2 ) 
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and that increases as magnitude of increases. For fixed' 
value of decreases as e increases and this effect is 

further enhanced for positive values of c^. (fig. 3.3). 

Effect of Ap is to increase V and in case of pressure drop 
(i.e. Ap < O) negative values of V are obtained (fig.3^4^3.5). 

Jtr 

Similar effects of Ap and c^ can be noticed frcwi fig. 3.6. 

Effect of 6 on V is shown in fig. 3.7 and it can be observed 

that an increase in 6 causes a decrease in V . 

P 

Effect O-f various paraneters on time average flux Q 

is graphically shown in figs. 3.8 - 3.10, In all the cases 

the effects of various parameters on Q are opposite to those 

on V^. It is noticed that 
P . . 

Q(ci,ei,£) = QC-Cj_/-e^/-e) 

It can be seen that Q increases with when e ^ 0, 
a reverse trend is, however, noticed for the case e > 0. 

A similar effect of c^ is also observed. For fixed values- 
of c^ and increase in 0 with increase in e is noticed. 

The effect of increasing Ap is to decrease Q. In the particular- 
case when e = 0, D decreases as magnitude of increase for 
Ap < 0 and an opposite trend is noticed for Ap > O (fig. 3.9). 

3.3 IRANSBDRT THROUGH A TUBE 

3.3.1 MATHEMATICAL MODEL (CASE 2) 

Let us consider the flow of an incompressible viscous 
flxiid in a cylindrical tube of radius a, where the central 
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core region (Fig, 3,12) of radius 6^ is moving with a 
velocity to be determined, in the negative axial direction'. 
The fluid flow is confined only to the region between the 


boundaries of core and the tube. Let a peristaltic wave 
of wavelength K and amplitude b' pass over the wall of the 
tube with a speed and let another peristaltic wave of 

amplitude b£ pass over the boundary of the core region. It 
is assumed that in the steady state, this peristaltic wave 
tries to align itself to the wave passing over the wall of 
the tube and hence acquires the same speed and wavelength. 

The shape of the wall of tube and the inner core region H£ 
are, therefore at the instant t' , given by 

H' (X',t') = a + b' sin - ct' + V£^t) (3 .31) 

H£(x',t') = 6' + b£ sin -^(x' - ct' + V^t' ) (3.32) 

—3 

The Reynolds number involved here is of order 10 
as in the section 3.2 ^ the equations of motion therefore 
are 


3X' 

dP^_ 

8R' 


+ 


+ 

^ (R' 

ay' 


R'c^R' 

dR* 

a^v' 

- ' ,.r -t* 

— (R' 

av' 

ax' ^ 


Sr' 


= 0 


= 0 


(3 -.3 3-) 
(3.34) 


5X' SR' R' 


(3.35) 
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FIG. 3.12 MOTION OF A MICRO ORGANISM THROUGH A TUBE 
{a)Without peristalsis (b) With peristalsis 
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The boundaries conditions are (as in section 3*2) 


U" = G'(X',t') = sin (X' - Ct' + V't') 

— X A D 


V' = 


6 H' 

at 


at R' = H' (X' , t' ) 

(3.36) 


and 


u' = -v; 


an' 

V' = --i 

at 


at R' = H'Cx'^tM 


(s-.s'T) 


The force equilibrium condition on the self propelling 
micro-organism is 


/ dS = AP' 

S 


(3.38) 


In wave frame moving with velocity o-V;^ in the positive 

hr 

axial direction, we define the following non-dimens ionalized 
variables 


X' - ct' + V^t 


£1 

a 


t'c 


P' 


X = - - 

U' _ c + Vl 
^ = - ■ — 

(g,Ci,Vp) = (g',c',vy/c. 


t-v-. _ 

,r-_ ,t - '^-,p - 


^ = Va ' ^h,h^,e,£^,6) = (H',H',b',b5_,6')/a 


The equations (3 .33)“ (3 *35) beccme 


_ 9E_ + p2 u 1 6 ( 9U ^ ^ Q 

_ -Sp g2rp2 + 1 8 .( 3y), , 0 

5r ^ ^ a 2 r 5r dr 


(3.39) 

(3 .40) 

(3.41) - 





9U ^ T_ 

as dr r 


= o 


( 3 . 42 ) 


where P = a/X. 


For long wavelength approximation P« 1. The equations 
C3.4O) and (3«4l)/ are therefore simplifed as, 

- r = 0 (3.43) 

||- = 0 ( 3 . 44 } 


Vie get the corresponding boundary conditions for u as 

(3 .45') 


u = g(x) + V -1 at r = hCx) 

ir 


u = -1 


at r = h^Cx) 


where 


g (x) = c^ sin 2 it x- 
h(x) = 1 + e sin 2 n x 
h^(x)= 6 + e 2 _ sin 2nx 

The force equilibrium condition (3.38) becomes 

1 - dh, (x) 2 

27t /[(rg|)^^. t h^(x)p]dx = Its Ap 

o i 

1 

and / p dx = O 
o 

Ap is the pressure rise over a wavelength- 
3 . 3.2 ANALYSIS 

The equation (3.43) can be solved, using boundary 
conditions (3.45) asr 


(3 .46) 


(3 . 47 ) 



u 


, 1 ^ r 2 ) (g(x)-i 

^ 4 ^ “ InCh/hp- - " - "li 


+V )ln(r/h, ) 
P ■ A - 
InCh/li^) 


(3 .48) 


as 


The flux in the wave frame, q{=Q^/i^ac) is obtained 


q = 2 / r u dr = i 


which can be rewritten as 


where 


with 


= |3 - 8 h3 


11 


11 


P P 


11 


Fi2(x) = h^-h^ 

2 ^1-^^ 

Fi3 (x) = h + 2rn’(h7h^T 
Integrating (3.49) over a wavelength we get 


Ap - 8Ij_^q - - ®^1^13 “ ®^14 

1 

♦ T 

12 


T 

^11 ~f p ‘M"'» f 
0-^11 


1 F^_ (x) 

A. = / p+%j--dx 


1 F (x) 1 

^13= ^ F,vur h4 =-f F'j-txJ 

o r j. ^ j--*- 


(3 .49) 


(3 .50) 


dx 


(3 .51) 


(3-52) 


Equation (3.5l) provides a linear relation in and q. 
Another relation is obtained from the force equilibrium 
condition (3.47) which, after substituting expressions for 
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“ (3-48), (3.49), beocme 


o = 4l2j^g - - Zl23Cj - 41^^ 


where 

^ o ^21 ^ 


1 FjjCx) 


• “T M ^ ^ 

' ^22 ■ { (xJ 


o ^21 


1 ^22^^^ 


X "22 ^ ^ (x) 

Ij3 = r t f.-Tx) 


and 5'2 ^(x) = (hj+h^)ln(h/h^) + (h^-h^) 

P,^(x) , 

F22(x) = 2 

Solving equations (3*5l) and (3.53) we get the 
expressions for V and q as 

Jr 




r^i3 

M 

oo 


q = 


^21 ^ ^22"=! ^ 


23 


M 


oo 


where M 


oo 




M 


11 


= 21 


21 


Mi 2 - 8(2I^3l23_ - ^23^11^ 
^13 ^ ^^^^ 14^21 “ ^ 24 ^ 11 ^ 


M 


21 


= I 


M 


22 


22 

® ^^13^22 - ^23^12^ 


^23 “ ® ^^ 14^22 “ ^^ 24 ^ 12 ^ 


(3.53) 

( 3 * 54 ) 

(3 , 55 ) 

(3.56) 

( 3 . 57 ) 

(3 .FS) 





The time — averaged flux C in stationary frame can 
be obtained by using following relation 

-P 2 e? 

Q = q + (1-Vp)(l .. ^2+ I - -i) (3.^9) 

3.33 RESULTS AND DISCUSSIONS 

The integrals given by expressions (3.52) and (3.54) 
are corruputed numerically using trapezoidal rule and the 
expressions for and Q (3.56, 3'.5"9) are ccmisited. It is 
observed that the qualitative behaviour of the results in 
this case is similar to those discussed in section 3.2.3 
corresponding to flow through a channel. 

In figs 3.13-3.18 , as a function of is 
plotted for various values of £ and Ap. is shown 
to decrease with e^, attains a minimum and increase again 
(a = + . 05 ) (fig. 3vl3 ) * For higher value of e (+ .2) 
this minima is not seen in the observed range of and 
monotonically decreases (a = .2) or increases (e = -.2) 

(fig. 3 . 14 ). Comparing these two figures (3.13,3.14) it 
can be observed that range of is increased as a increases. 
The effect of increasing c^ is to increase the magnitude of 
V and that of Ap is to decrease it. The maximum (or minimum 
negative) value of V is obtained wben 0 arxi Ap y 0» Effect 

Jtr 

of 6 is also similar as observed for channel case for 

( let 5. *05) (fig- 3.18) ♦ 
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FIG 3-19 EFFECT OF C] AND € ON Q (TUBE ) ( 8 =0-2 , ^p= 0 ) 
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FtG. 3.21 EFFECT OF fip AND Ci ON 0 (TUBE )( 6 = 0.2 ,€ = 0 ) 
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Similarly the effect of various parameters on 0 
can be seen from fig (3 .19.)-(3 .22 ) which is similar to 
those discussed in case of flow through channel. 

3 *4 application TO OVUM-TRANSPORT I FERTILITY CONTROL 

The ovum is the fanale germinal cell and is one of 
the largest cells in human body (diameter of ovum ' 150-200 Jim). 
At the time of ovulation the ovum is extruded from the ovary 
into the ovarian end of the fallopian tube and is immediately 
sucked in by the transverse; motion of smooth muscles of the 
finger-like projections existing there and by the cilia 
beating in waves. As the life span of the unfertilized 
ovum is very short (approx. 18 hours) it must travel fast 
enough to reach the other end of the oviduct to get fertilized. 

It is thus evident that ovum transport along the oviduct is 
a vital link in the process of reproduction. In recent years 
there has been enough interest in the study of biochanical 
and physiological aspects of this problem. A review of such 
investigations has been presented by Pauerstein and Eddy (l979). 
As suggested by experimental data, cilia are necessary for 
ovum transport and there appears to be a correlation between 
the number of ciliated cells on the fimbria and fertility in 
humans [^Brosens and Vasquez (1976)/ Pauerstein and Edc^ (1979)]. 
The muscular activity in the isthmus of the oviduct is also 
observed to be related to ovum transport [Hodgoson et al.(l977)]. 



\ \ 



3 23 EFFECT OF C, AND Ap ON VARIATION OF Vp WITH 
fR-0.2 ) (OVUM TRANSPORT THROUGH A TUBE) 
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Recently Blake et al. (1983) studied the effect of 
muscular and ciliary activity along the wall of the oviduct 
by adding a force distribution term in the equation of motion- 
However, these motions can also be accounted for by assuming 
a combination of transverse and longitudinal motion of the 
walls as discussed in the above models. Since there is no 
evidence of ovnm having its own motility, a model of the 
ovum transport through oviduct, accounting for the muscular 
and ciliary activity of the wall and the pressure difference 
across the oviduct, can be obtained by setting = O in the 
model discussed in section 3.3. 

In this case the propulsion velocity has been 
computed using expressions (3.56) and (3.58) by putting 

= 0. These values of are plotted as function of £- in 
fig. 3.23 and the effects of Ap and c^ are graphically sho^'. 

It is observed that becomes more negative as the magnitude 
of e is increased (i.e« the ovum is carried over by the 
peristaltic waves). It is also noted that increases as 
Ap and c^^ decrease. 

Since the ovum is carried along by the peristaltic waves, 
it is observed that its propulsion speed can be reduced if 
the amplitude of these waves are either minimised or are applied 
in the direction away from the uterus. It can therefore be 
suggested that for fertility control the ciliary and muscular 



Ill 


activity should be as minimum as possible which could.be 
induced by using some biochemical means such as hamone^. 

3.4 CONCLUSION : 

Models for swimming of a micro — organism through a 
channel and a tube have been proposed , taking into acco\int 
the effects of peristaltic motion of the wall and the cilia 
motion* This model has been applied to ovum transport 
and it has been shown that the velocity of the ovum towards 
the uterus can be reduced by a reduction in the cffnplitude 
of peristaltic and longitudinal waves- This reduction in the 
propagation speed of ovum is important frcsn the point of 
fertility control# due to limited life span of unfertilized 
ovum. Thus if the contraction and relaxation of anooth 
muscles of the fallopian tube, causing peristaltic wavesr 
and ciliary motion can be controled by suitable combination 
of harmones, the chances of fertilization can be considerably 
reduced. This aspect should therefore be explored further^ 
both theoretically and experimentally. 
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CHAPTER IV 


EFFECTS OF VISCOSITY VARIATIOH OP A 3I0FLUID ON THE 
MICRO-ORGANISMS UNDER PERISTALSIS 

5.1 INTRODUCTION 

^The study of self propelling micro-organisms was 
initiated by Taylor (1951) and has since than attracted attention 
of many researchers [Hancoch (1953), Shack and Lardner (1974), 
Pironneau (1975), Shukla et al. (1978)]. In chapters II and III 
of this thesis this probl^i has been studied, taking into 
account the peristaltic and longitudinal motion of the walls 
of the transporting ducts which may arise due to muscular and 
ciliary activities therein* These models have also been applied 
to study the motion of spermatozoa and ovum* effect of 

finite thickness of the micro-organism on its motion has also 
been studied in chapter III. 

In these studies the effect of viscosity variation of 
the fluid in which the micro-organisms swim has not been 
taken into account. However it has been suggested in the 
litrature that the peristaltic transport may be involved in 
sperm transport through vas deferens [Guha et al* (1975“), Gupta 
and Seshadri (1976) Srivastava et al. (1983)] and in tJ^ 
cervical canal [Smelser et al. (1974)] - all being ti^ cases 
where the viscosity of the fluid near the wall is different 
from its viscosity in the center of the (Sict* Xn view of this, 
a model has been presented by Shukla et al- (1980) to stiMy 
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the effects of- viscosity variation of the fluid on the 
mechanism of peristal-tic -transport through a pipe and a 
channel* It niay be noted here that in these studies "the 
interaction of peristaltic waves generated by the micro—organiarn 
on its body/ to s-wrim across the duct/ -with the viscosi-ty vstria-tion 
of the bio—fluid has not been s-tudied. It is known that the 
cervical mucus resembles like a close mesh during the luteal 
phase and is difficult to be penetrated by the Spermatozoa 
[Odeblad (1959), (1962), (1968)], Ifowever, during the mid-period 
of the menstrual cycle the micelles are aligned near th« wall 
[Shukla et al, (1978)] and a peripheral layer having different 
viscosity is formed near the wall of the cervical canal. 
Investigations has been made to utilise this behaviour of 
mucus as a means of contraception [ James and Marriott (1983)]- 
In view of the above, mathematical models have bean 
presented, in this chapter to study the transport of a 
micro-organism in circular and rectangular ducts- under the 
effect of peristaltic and longitudinal wave motion of -the 
walls and the viscosity variation of the bio-fluid along 
the cross-section. The equations of mo-tions and the force 
equilibrium conditions applied in chapter III are iiK>dified 
to take into account -the visvosity varia-tion. As in 
chapter III the solutions have been obtained for, this- gmeral 
case, A particular case when viscosity is a step fxonction 
has been discussed. It has been shown that, apart from 
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the results obtained in chapter III, the propulsive velocity 

of the miciro— or'go.nisru decreases as the viscosity of the 
peripheral layer increases. 

4.2 TRANSPORT THROUGH A CHANNEL 

4.2.1 MATHEMATICAL MODEL 

Consider the propulsion of a flexible slab 
( micro-organism ) of infinite width and finite thickness 
25 ^ in the middle of a channel of thickness 2a. The space 
between boundaries of the slab'and the channel is filled 
with a Newtonian fluid with varying viscosity. The walls 
of the channel are vindergoing both peristaltic and longitudinal 
wave motion. It is further assumed that the slab is being 
propelled at a constant speed in the negative X' -direction, 
in a stationary frame (X'.,Y'), by propagating peristaltic 
waves along it surfaces. If these waves travel with the 
speed c-V^ in the stationary frame then in this frame the 
shapes of the walls of the channel (Y* = + K' ) and the 
slab (y' = +Hp at any instant t' are given by (see fig- 4'.l) 

H' (X' ft*) == a + b' sin ^ (X* -ct' tV^t* ) ( 4 .I) 

Hj^CxSt') a 6" + b^sin ^(X'-ct'tV^tM (4.2) 

where b' and b' are the amplitudes of waves on boundaries 
of the channel and the slab respectively. 
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FIG.4.1 FLOW OF FLUID IN CHANNEL WITH PERIPHERAL LAYER 

(a ) Three dimensionat without 

(b) Two dimensional with peristalsis 



fhe equations governing the flow of the fluid 
through the gap Y'< H'), under long wavelength approxlmat- 

ion and for low Reynolds number condiUons, can be written 
as 


O 

0 


+ -8-(n' -W'.) 

3X' 3Y" 3Y^ 


3^ 

3Y' 


(4-3) 

(4.4) 


a 

ax' 



(4>5-) 


where U' and V' are the velocity components in X^ and Y' 
directions respectively, P* is the pressure and = 1 ^* (Y' ) 
is the viscosity of the fluid which varies along the cross- 
section of the channel. 


i.e. 


The boundary conditions are sane as in section 3.2, 


U' = G' (X'-ct'+V't') = c' sin ~.(x'«ct'+V't) 

p - 1 A p 


U' = 


at Y' = + H' (X' , t' ) 


at Y' = + (X' , t' ) 


(4.6) 


Force equilibrium condition on the surface of the 
micro-organism is 


(4.7) 


/ T dS = 6' 

S 

where x is the resultant of the force acting 
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surface S. of the micro-organisir. and Ap is the pressure 
rise over a wavelength. 


Applying following transformation 


X = (X--ct'+V^t)/X , y = Y'/a, t = ctVA 

u = (u'-ctV^)/c , V = W'/ac, tgrc^fV^) = (g',c£,v^)/c 

p = p'aV^ir 4 ' ^e»ei/6) = (b%b£,6M/a 

where = 4 ' (0),the eqn, of motion {4.3)-(4.5) can be 
written in dimensionless form of v/ave fraine moving with 
a velocity c-V*, as 


0 - „ + -3 . (u SLli) 

3x ^ 3y ^ ay^ 

0 = - e 

37 


+ av ^ 0 

3x ^ 37 


The boundary conditions 6) are transformed as 


u 


g (x) + V -1 at y = h(x) 


= -1 


at y = hj^(x) 


where 


h(x) = 1 + e sin 2 7 Tx 
h^(x)= 6 + sin 2 n x 
gtx) = sin 2 ti x 


The force equilibrium condition, in this case becannes 

= 6 Ap 


1 

/ 

o 


i7>y=h,(x) ^ 


dhj^ (x) 


dx 


(4.8) 

(4.9) 

(4.10) 
{4'.ll) 

(4.12) 

(4-13) 


(4.14) 
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4-. 2-. 2 ANALYSIS 

The equation (4.9) can be solved, using boundary 
conditions (4.12), as 


where 


= ( fftj, 


I,(y) 


1 
h 


4 


4)- dy 


“1 

I, J axlr 


(4 -.15) 


(4.16) 


The flux in the moving frame is obtained as 

h l5(h) I,,(h) 

q = / u dy = -(h-h^)-(|J)[l 3 (h)- f (y] + (g+V )[h- 

1 11 

(4.17) 


which gives as 

3 X 


) = 




2"* f {q+(h-h^)} l3^(h) + (g+Vp){l2(h)-hl3^(h }] 


3 

\ ji * 1 

2 

(4.18) 

Integrating 

(4.18) 

we get 



-Ap = 

^ll'3 + ^12^p ^13^1 ^ ^14 

(4.19) 

where 

^11 

1 I. (h) 

f [■-, 2- 

o I 2 _(h)l 2 (h)-l 2 ^h) 
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12 


13 


=. / . ,^2(h)-hli(h) 

° IiCh)l3(h)-l2(h) 
1 ^ I^(h)-hl(h) 

= I r -^ 2 , _ 


] dx 


2 ^ ]sin2Tr x dx 


(4.20) 


1 ^ (h-h, ) It (h) 

^14 = ^ I - _ ]dx 

o ^iChJia Ch)-.i2(h) 

Using expressions (4.15 ) and (4.17) for u and p, the 
force equilibrium condition becomes 


where 




■21 


= ^ [ 


I^Ch) 


*-]dx 

o Ij_(h)l 2 (h) - 12 (h) 


1 I, (h) - hl.(h) 

? = / ['■ ■ > . . . - -jdx 

o Ij^(h)l 2 (h)-l|(h) 


1 13 (h) - hl 2 (h) 


] sin2 7i X dx 


(4.21) 


(4.22) 


^ ' -- 0 

o I^(h)l3(h)-l2(h) 

^ (h-hT)l.(h) 

o ‘‘ I 3 ^(h)l 3 Ch)-l 2 (h) ■ 

The equations (4.19) and ( 4 . 21 ) are linear in q and 
and can be solved, therefore, as 

MiiAp + ^^2*^1 ^13 


= 


M 


00 


M. 


( 4 . 23 ) 


00 



where 


M 

oo 

r: 

F F 

22 11 

- P P 
12^21 

^11 

= 

F 

21 


^12 

= 

h3‘’21 

- F P 

23 11 

M 

13 

= 


- P P 

24 11 

M 

21 


F 

22 


M 

22 

= 

F F 

13 22 

_ p p 

23 12 

M 

23 

= 

F F 

14 22 

- P P 

24 12 


The time averaged flux Q in statiaiary can be 
obtained from the following relation 


(4.24) 


Q = q t (1-V^) (1-0 (a. 2^) 

^ p - 

Eqns . (4.23 )- (4. 25) determing the propulsion velocity 
Vp and the time averaged flux in the fixed coordinate 
system for any viscosity function |i. (y) . In the following 
section effect of a peripheral layer is considered as a 
particular case when 4 (y) is a step fmction. 

4.2*3 EFFECT OF PERIPHERAL LAYER 

To Study the effect of peripheral layer the viscosity 
function has been assumed in following form 


= il / “h < y < h 

jd (y) =r 

= 1 / y.lcch 


(4.26) 


where (l-a)h is the thickness of the peripheral layer 


near the wall of the channel 



VJxth this viscosity function^ the integrals 
I^(h), l 2 (h)^ ^ 3 ^^^ can be analytically evaluated as follows 
Ij_(h) = ^(h-ah) + (ah-h, ) 

12 (h) = -^(h^-A^) + i{a\^-h?) (4.27 

2|ix ’^1 

(h) = “^(h^-a^h^) - 1 - i(a^h^-hT) 

"^34 2 1 

It can be seen from (4.27) that this model reduces 

to the model discussed in section 3.2 of chanter XII^ when 

there is no peripheral layer i.e., -a = 1 or 4 = 1. 

4.2.4. RESULTS AND DISCUSSIONS 

Assuming 4 (y) as a step function (eqn. (4.26))# to 

account for the presence of peripheral layer# and 0 have 

been calculated from expressions (4.21), (4.23 )-(4. 25)# (4.27) . 

These results are graphically shown in fig. 4 . 2 - 4 . 7. It 

can be observed that effects of Cj^ and e on (fig. 4 . 2 ) 

and Q (fig. 4.5’) are similar to those obtained in section 3.2 

of chapter III, i.e. V decreases and Q increases as c and 

CjL increase. Effect of 4 and a on and Q can be sti:uiied 

from figs. 4.3 and 4.6 respectively* It is noted here that 

magnitudes of V and Q decrease as oc and 4 increase. Effect 

P 

of Ap is also similar as observed in chapter III i'.e. as Ap 
increases from negative to positive values. beccanes more 

negative (fig. 4-.4') and S becones more positive (fig. 4.7). 
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c.r / 9 PFFFrr OF Cl AND Z ON VARIATION OF Vp WITH E, 
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FIG. 4.4 


EFFECT OF u AND ip 01«l VARIATION OF Vp WITH C, 
X / r- A o r --2 £ =— 2 , «C = -90) 

(CH ANNEL )( 6=0*2, C| - - a 
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4.3 TRANS EORT THROUGH A TUBE 

4.3.1 MATHEMATICAL MODEL 

Consider the flow of an inccmpressible Newtonian 
rluid/ whose viscosity varies along the cross-section of 
the tube/ in a tube of radius a^ where the central core 
region (fig. 4 .8) of radius 5 ^ is moving with a velocity 
Vp in the negative axial direction. The fluid flow is 
confined only to the region between the boundaries of the 
core and the tube. Let a peristaltic wave of wavelength K 
and amplitude b' pass over the wall of the tube with speed 
c-V' . It is also assumed that a peristaltic wave of senje 
speed and wavelength and amplitude b' passes over the 
boundary of the cylinderical core. Then the shapes of the 
boundaries of the tube and the core , at instant t' 
are given by 

H'(X',t') = a + b'sin ~(X' -ct'+V' tO (4.28) 

A P 

H5_(x',t') =6' + b' sin ^(X^-ct'+V^t') U.29) 

Foxr low Reynolds number the equations governing 
the flow, under long wavelength approximation are given by 


a p' 

t i- 

J - 

OXR' 3 "',.) 

(4.30) 

"a X' 
ap'. 


aR' 

aR' 

(4.31) 


a R' 

^ 1 _ .9 . ) = 0 

ax' R' aR' 


(4.32) 
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Cl nm IN tube 

)W Oi- etnists peristalsis 

I Without p^nstalsts 
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Considering the longitudinal motility of the wall 
of the tube the boundary conditions are 

U' = G£(X'-ct'-!-V^t') = c^sin -ct'+V' t' ) 

at R' = H' (X^t') (4.33) 

U' = -V^ at R' = H£(x',tM 

Transforming all the quantities to wave frame 
(moving with velocity c-V' in the positive axial direction) 

Jr 

and non-dimensionalising them with following relations 
X = (x' -ct'+V^t' )/\ , r = R'/a # t = t'c/^ 
u = (u^-c~V^) , 
v = V' ?y'ca, p = p* 


(e,e,6) = (b',b^,6M/a 

(4.31) 

(where is some characteristic viscosity/ 

the equations 

(4 .30) -(4. 32) can be written as 


0 = - 2lE + 1 ^ (4r^) 

3x r ar 3r 

(4.35) 

0 

tl 

1 

(4.36) 

^ + aji + Y = 0 

9x a r r 

(4.37) 


and the boundary conditions (4.33) are transformed as 

u = g + Vp-l ; r = h(K) , 

= -1 ; r = hj^(x) 
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where 


g(x) = sin 2n x 
h(x) = 1 + e sin 2n x 
hj^(x)= 6 H- sin 2n x 

The force equilibrium condition, in this case is 

1 r / u X dh, (x) - 

27T / [ + h 3 _(x) -gj-' P]dx =n6^Ap 

o 1 


(4.39) 


/ p dx = 0 
o 

4.3.2 ANALYSIS 

The equation (4.35) can be solved, with bomidary 
conditions (4.3 8), for velocity u as 


(4.40) 


I. (h) 


I^tr) 


T J-T Vll/ * 

u . -1 - A-h, I^Cr)] + (9+Vp) Ajj 


where we define 


(r) 


= / 


irCr) 


dr 


(4-41) 


- / 


dr 

idTrT 


(4.42) 


r ^3 


^ IJ.tr) 


dr 


h. 


The dimensionless flux in the wave frane jq - O'/ ^ a cJ 
is obtained as 

h 

f. = 2 1 ru dr 

h. 
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which using eqn, (4.41) gives 
"J = - 13 (h)- ftl-hVi-th^-hJ) + (g+Vp)[h=- 

2 

and therefore x^re get 


_ i(l£) 

2 3 X 

1 ^9 

= ''V 

( 4 . 43 ) 


ihh) 


where 

°l = "3 - TJ(hT 

(4.44) 


P I, (h) 

^2 = ^ - i^Th) 

Integrating (4.43) we obtain the follovang relation 


in q and V 


- 2 “ ^^11 ” ^ P ^12 " ^ 1^13 ^14 

1 


(4. 45 ) 


with 


11 


. dx 

o ' "^7 




12 


13 


14 


t 


( 4 * 46 ) 


= ^ [ 


G_sln 2 71 X 


- ]dx 


1 h^-h^' 


The force equilibrium conditions (4.40) provide 
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where 


I, (h) 

1 Gt| — G I- (ll) 

; r Idx 

n n. I _ (h) 


0 Gj^l2(b) 

^ G,-G„I, (h) 

F23 = / [ -S-r^Th-r^i" 

1 h^-h^ 

^24 = ; 


The equations (4.4S) and (4.47) can be solved 

%.jr 




^11 ^12'^! '*’- 


^21 ^2_2fl ,^23 . 

q = ■■ '“ '" m' V 


where 


^00 = ^21^12 ^ ^22^31 


" 2^21 


Mt2 = --21^13 " ^23^11 


^13 = ^2/14 ' ^2/11 

1 F 

^21 = ■ 2 ^22 

M22 = ^22^13 - ^23^12 

M23 = “^22^14 " ^24^12 


(4.48) 

(4.49) 

(4.50) 

(4.51) 
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The time average flux in stationary frame can be 
obtained from the following expression 


2 £ ^ 

Q = q + ( 1 -Vp)(l - 52 + ^ 

4 . 3.3 EFFECT OF PERIPHERAL LAYER *. 


(4.52) 


To Study the effect of peripheral layer, the 
follovri-ng viscosity function has been considered' 

M-Cr) = 4 t ah< r< h 
= 1 f 0 _<r_<ah 


(4.53) 


With this viscosity function the integrals 
Ij_(h), 12 (h), 13 (h) are evaluated as follows 


I, (h) = l(a^h^-h?) + — Ch^-aV) 

1 2 1 

12(h) = 


(4.54) 


I- (h) = l(aV-hf ) + (h^-a'^h^) 

^ 411 


It is evident from eqn. (4.54) that the proposal 
model reduces to the model discussed in section 3 .3 of 
chapter HI, when there is no peripheral layer’ i.e. 
a = 1 or /I = 1 . 


4.3.4 RESULTS AND DISCUSSION 

The expressions of end Q have been oasputea frcm 
equations ( 4 . 49 )-( 4 . 5 - 4 ) ana are gralMcally showi in 
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FiG.4.12 EFFECT OF AND E ON Q (8=0-2,iip= 0 

jd s 0.75 r ^ =0*90 ) (^^8E ) 
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RG.4.,4 £FFECT Of and Ap ON RATION Of Q WITH C, 

(6=0.2,C,=-.2,e=-Ob,^ = -90HTLBtJ 



figs. 4.9 - 4.14. Comparing the results with those 
obtained in section 3.3 of chapter III, it is noticed 
that for fixed viscosity, effects of e and c, on V and 

1 P 

Q remain same. The effects of |I and a are qualitatively 
similar as observed in section 4.2. i.e. magnitude of 
Vp and Q decrease as h decreases. 

4.4 CONCLUSION *. 

Mathematical models to study the effect, of peristaltic 
and longitudinal motion of the walls of the transporting 
duct and the viscosity variation across the cross-section 
on the swimming of a micro-organism have been proposed. 

It has been observed thab the results obtained for the 
transport in channel are qualitatively similar to those 
of the tube. Apart frcati the results observed in chapter III, 
it has also been shown that magnitude of propagation velocity 
of the micro-organism decreases as the viscosity of the 
peripheral layer increases. 
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CHAPTER V 


EFFECT OF MULTIPLE CONSTRICTIONS ON PERISTALTIC TRANSPORT 
OF A FLUID THROUGH A TUBE OF NON-UNIFC®M GEOMETRY 

5 . 1 INTRODUCTION 

is a mechajiijOT ior the transport of a 
fluid achieved by the passage of progressive waves of 
contraction and expansion along the walls of a distensible 
tube. Many physiological fluids in the human body are 
known to be transported by peristaltic action of the walls, 
such as urine in ureter and chyme in intestines. It is 
also speculated that peristalsis may be associatei with 
the vasomotion of small blood vessels [Fung and Yih (1968)] 
Severals authors have contributed to the study of 
peristaltic transport in both physiological arxi mechanical 
situations Burns and Parkes (1967), Fung and Yih (1968), 
Barton and Raynor (1968), Shapiro et al. (1969), and 

Ostrach (1970), Manton (1975), Kaimal (1978), Shukla et al, 
(1978)]. In particular- Rung and Yih (1968) investigated^ 
this problem for arbitrary wavelengths and small amplitudes 
Later, Shapiro et al* (1969) considered peristaltic 
pumping in a tube under long wavelength approximatiMi 
and for very small Reynolds numbers. It has been shown 
in these studies that the average flux over a time period, 
is increased by increasing the amplltu(te of the peristaltic 



wave. Li (1970) analysed peristaltic transport in 
circular cylindrical tube, assuming stream function 
as a power series, and canpared the flow through tube with 
the flow in two dimoisional channel. A detailed analysis 
of peristaltic pumping using finite element method has 
also been presented by Tong and Vawter (1973). Manton (l97s) 
has studied this problem by perturbation techniques, using 
long wavelength approximation .The effect of non— Newtonian 
behaviour of the fluid under peristaltic motion has also 
been considered by many investigaters [Raju ^d Devanathan 
(1972), Acharya (1982), Shukla and Gupta (1982)]. 

It must be observed that in these studies the solutions 
have been obtained in the wave frame, vdiere the flux is 
steady because of the uniform geometry of the transporting 
ducts. However, it is noted that cross section of various 
ducts in human body gradually varies along their length and 
therefore a steady flow through them can not be obtained 
even in the wave frame. Guha et al. (1975) studied the 
effect of peristalsis on the flow through a duct of linearly 
varying cross — section with particular reference to 
spermatic fluid transport in vas deferens. Later Gupta 
and Seshadri (1976) and Srivastava et al. (l983) also 
studied this problem and obtained the solution in the 
stationary frame* In these studies the value of pressure 
difference across a wavelength has been calculate! for ar 
given rate of flow and periodic nature of the pressure 
difference has been shown. 
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Further, it may be pointed out here that the peristaltic 
transport of fluid in constricted tubes has not been studied 
so far, though such situation may arise in intestines due 
to the development of polyp. The dynamical interaction 
between the passage of the progressive peristaltic waves 
and an abnormal growth in the lumen of the transporting 
duct may be very complex' and its study may lead bo a better 
understanding of such pathological situations* In view- 
of this, in this chapter a mathanatical model has been 
proposed to study the .peristaltic transport of a Newtonian 
fluid through a circular cylindrical tube of varying 
cross — section with two mild local constrictions. Following 
Gupta and Seshadri (1976), solutions are obtained in the 
fixed laboratory frame under long wavelength approximation 
and the effect of various parameters on the pressure 
difference and the wall' shear stress have been studied. 

5.2 MATHEMATICAL ANALYSIS: 

Let us consider the flow of a Newtonian viscous- 
fluid through a circular tube of non-uniform cross-section 
and having two constrictions along its length. The geometry 
of the tube, as seen in the stationary frame of reference, 
is shown in fig. 5,l-> is assumed that sinusoidal 

travelling waves are imposed on the flexible bourxlary of 
the tube. The radius of the tube, at a distance z' from 




FVg.S'l Geometry of the tube with peristalsis, multiple 
stenosis 
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origin and at instant t' is given by 

R (z / t ) = RgCz') + b' sin ^(z'-ct') (5vl 

where b' , X and c are the amplitude^ wavelength and 
speed of the jJeristaltic wave and (z^ ) is the variable 
radius of the lumen of the tube when there is no peristalsis. 

The constrictions are assumed in series^ mild and 
are symmetric about axis and therefore the lumoi of the tube 
can be assximed to be given by 


I 




0 < z' < d' 


^o “ ~Cz'-d'- -|)];d£< z' £ d£<L' 


R' Cz' ) 
s 




;d£+L^£ z' < - 2 


L' 

2 


a; 


> - ^[1 + oos ^Cz'-B' )] ;b'- <B^ 

6 ' K 

■ R’*'(z' ) - -f [1 + OTS ^(z'-B')]; B' < z'< + -f 

■ ^2 

i K 


(5.2) 

As shown in fig. (5.l)/ the constrictions are in the regions 
d£< 2 '< d'^+L' and B£-L'/2< z' < B5_+L'/2 respecUvely.R' is the 
constant radius of the unstenosed part of the tube in the region 


*^11 2'< and R''''^(z') is gradually varying radius in the 
region B£< z'< B' . and 6^ are the lengths and the maximum 

thickness of the two mild stenoses (the suffixes 1 and 2' 
indicate the first and , second stenosis respecUvely) and 
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are such that they follow the restrictions for mild stenosis 
[Young (1968)] i.e'. 

r 6d« min (R',R^ J 
i out 

(i = 1/2) (5.3) 

<< li' 

iL- 1 1 

where 

The flow is considered to be inertia free, laminar 
and axi ally-symmetric. Following Young (1968) and Shukla 
et al» (1980) the basic equations governing the flow are 



(5.4) 

a z' r' dr' dr' 


9wi + 1_ L..(r'u') = 0 

is.s) 

a a' r' ar' 



where M- is the viscosity of the fluid, p' is the 
pressure and u^ and w^ are the velocity cwnponents of the 
fluid in the radial and axial directions respectively. 


The boundary conditions are 

w' = 0 at r" = R' (z^t') 


(s.e). 


ar' 


0 


at r' = 0 


Using the following non-dimensional quantities 


-v/D 2 


t == Ct'/K w = w^/c, P = 

(z/%/Ljl,L2,Bj_) = (z',d5_,b',L',Bj^)/B' 

(r,R,R^,R*£,6i)= (r',R',R',R*'-b''«P/®; 


(5.7) 
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and solving equations (5.5) and (s'.e) we have 

w = - (r2_^2j 


(5.8) 


The Instantaneous volumetric flow rate, in dimensionless 
form can be written as 

ri/’ ^ ^ t)]^ 

Q(z,t) =27 r w dr = - - - .,L4J9 

O ° O 2 

which gives 


,12 
9 2 




[E(s,t)]'' 


(^^.9) 


The eqn- (5.9) can be integrated with respect to z, 
over a wavelength to give instantaneous pressure rise over 
a wave length Ap(t) as 


Ap(t) = .-8 




(5.10) 


The pressure rise over a wavelength has been calculated 
for a given form of Q(z, t) as the inverse problem is difficult 
due to non-availability of analytical solution of the integral 
involved. 


It is evident from eqn. (5 .10) that Ap(t) will always 
be negative for a constant q(z, t) implying no pumping action. 
Therefore Q(zlt) has been assumed to be a function of (z-t) 
of the following form [Gupta and Seshadri (1976) J 


1 „2 


q(z, t) =□+[,£■ (z#t)-Rg(z/ 1)] ” 2 ^ 


( 5 -. 11 ) 



where Q is the constant time 


]5o 


mean flow over one period of 
the wave. It should be noted here that the flux in wave 
frame gCz, t), as defined by Jaffrin and Shapiro (l97l), 
will be a function of z in this case. 

The equation (5.10), therefore, beccmes 


Ap(t) = -8(q- 




1 

/ 


dz 


o [R(z, t)]' 


+ 8 


1 r2(2) 

/ '• ' dz 

o [R(z, t)] 


-8 / 


dz 


(5.12) 


o [R(z,t)]^ 

which gives the pressure rise Ap (t) in terms of the average 
volumetric flow Q. 

The instantaneous shearing stress t^(z, t) on the 
wall of the tube is given by [Manton (l97l)] 


(z* t) 




D(z,t)[r -f-J] 

L ar-i^ ^ 


(5.13) 


where 


D(2,t) = 


p2 ^dR^ 2 


ltp2(f)2 


(5.14) 


and P = R^/B' . (S'- IB-) 

Integrating (5.15) over a time period and using (5.14), 
(5.15) alongwith the form of Q(z,t) from (5.11) we get the 
time-averaged shear stress t(z) on the wall of the tube 
as 

f \ f- 1 2 d 2/ dt - 4 / D(z,t)dt (5.16) 

t(z) «-4(q- -^ s ^ o tR(z,t)f o 
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5.3 RESULTS AND DISCUSSIONS ; 

To explicity see the effects of various parameters 
on the pressure rise and the shear stress, the expressions 
for Ap and t^(z, t) on the second constriction ((5.12), {5.16)) 
have been numerically calculated on DEC-1090, The integrals 
involved , are computed using trapezoidal rule which provided 
most accurste results (accuracy upto 10 ^) in ccxiiparison 
to some other quadrature formulae. To study the effect 
of varying cross-section, R’^(z) has been taken as 

r’^(z) =s exp.[K(2-Bj^) (S'.IT) 

where the parameter K determines the convergence/divergence 
of the tube, the negative, zero and positive values of K 
correspond to converging uniform and diverging tube 
respectively. 

The following values of the parameters regarding the 
constrictions have been chosen in the calculation for Ap 
and T 

d, =0,2, L, =0,2, 1^2 = 0,2, = 0,8 

dm ^ 

and P = 0 , 1 . 

Values of Ap given by eqn. ( 5 , 12 ) are plotted as a function 
of time t for various values of Q, K and 5^,62 in (fig.5.2- 
5.4) . For fixed values of 6 ^^ and 63 / effect of Q 
and K is shown for e =* 0,4 in fig* 5'.2 and e - 0*5 in 





2 - 60 j- 

2-5oL /\ 

0 ' 



Fig. 5*4 Effect of Si and 5i on variation of 
with t (€ ='4 , Q = 2 , K = 0 ) 



K and e, the effect 


S'* 3/ while for fix«i values- of Q , 
of and is shown in fig, ( 5 , 4 ), 

For fixed values of Q, K , e , 5 ^ and 63 / oscillatory 
behaviour of t with respect to t can be observed (fig, 5 . 2 , 5 . 3 ) 
As t increases Ap goes down, attains a minimum at about 
1/3 rd of the time period of peristaltic wave ; it attains- 
a local maximum at t, slightly greater than half time period 
and the second local minimum is obtained about 3 / 4 th of 
time period. It must be noted that this behaviour of two-plus 
one extrema i.e. one maximum between two local minima is 
due to presence of two constrictions and is absent in the 
results obtained by Gupta and Seshadri (1976), Effect of 
increasing the flux is to decrease pressure rise Ap(0 = 0 
and 0.2) and it beccmes negative for 0 = 0,6 (see.fig* 5*2) . 
Also comparing figs. (S'* 2) and (5'. 3') it can be noted that 
Ap increases with increase in e and decrease in 5. Effect 
of increasing K is to decrease the magnitude of Ap. Comparing 
with the results obtained by Gupta and Seshadri ilSlS) it 
can be observed' that (-)ve values of Ap can be obtained 
at a lesser value of 0 in the proposed model and this effect 
may be due to the presence of constrictions. 

Fig, 5.4 shows the effect of 62 _ 63 

absence of any constriction = 0 / ^2 ~ presstire 

attains a maximum and then reaches its minimum value. For 






a single constriction (63^=0,62=.!) the local maxinium 
disappears from the given range. One maximim and two 
minima are observed when both constrictions are present 

(63^=-1/52=*1^/ the maximum being much higher than any 
Other value* 

Effects of Q and K are shown, in figs. 5". 5 (e = , 4 ) 
and fig. 5 . 6 (e = ,5), on ^2/ the stress on the secorsi 
constriction ( I z— L2). The maximum value of T2 is 
observed at the point of maximum thiclcness of the constriction. 
The magnitude of T2 decreases as Q increases and increases 
with e and ^ 2 ' 

5 . 4 ' CONCLUSION : 

A model has been presented to study the peristaltic 
transport of a Newtonian fluid through a circular tube 
of varying cross section and having two constrictions. 
Following Gupta and Seshadri ( 1976 ), the solutions have 
been obtained in the laboratory frame and it has been shown 
that the pressure rise over a wavelength and the shear 
stress on constrictions increase as the amplitude of the 
peristaltic wave and the thictaiess of the constrictions 
are increased. This res^t can be applied to the case o£ 
intestines having abnormal constrictions. 
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CHAPTER VI 

FLOW OF NON- NEWTONIAN FLUID THROUGH A TUBE WITH 
MULTIPLE STENOSES 

6.1 INTRODUCTION I 

It is well known that stenosis, which refers to abnormal 
growth in the lumen of artery, is a result of esi ar'terial 
desease. It is caused mainly due to intravascular atheroscl- 
erotic plaques which develop at the arterial wall, thus 
occluding the lumen of_the artery partly or fully. Although 
the etiology of the initiption of the stenosis is not clearly 
understood, it has been realized that hydrodynamic factors 
such as wall stress, pressxire distribution and separation 
phenomenon play an Important role in the developnent and 
progression of this desease [ Rodkiwitcz (1977) and^ Stehbens 
(1982)]. Hence considerable attention has beoi paid to the 
theoretical and experimental studies of blood flow in stenotic 
region under different conditions [ Yoxing (1968), Caro et al. 
(l97l), Nerem (1974), Deshpande et.al. (1976) , Me Donald (1979), 
Shukla et. al. (I980a), Padmanabhan and Devanathan (1981) and 
Talukder et al.(l976)]. Most of these analytical studies 
have investigated the flow characteristics of blood in arteries 
with mild/non-critical stenosis^ by considering blood as a 
Newtonian Fluid [ Young (1979)]. However it may be points out- 
that blood being a suspension of cells in plasna, has=be^ 
found to behave as a non— Newtonian fluid at low shears rates* 



in tubes of small diameters. In view of this some investigators 
have considered non— Newtonian behaviour of blood [Shukla et al. 
(1980b and c), Williams and Javadpour (1980), Devanathan 
Parvathamma (1983)]. Shukla et al. (1980b) studied the flow 
characteristics for mild stenosis by considering a non-Ne%rt 3 onian 
power law fluid model for blood. They have also considered the 
effect of peripheral layer on the flow of blood to take into 
account the clustering of cells in the coatral zone ^ Shukla 
et al. (1980c)]. 

In above mentioned studies only the effect of single 
stenosis has been considered and the tube waa taken to be of 
xiniform cross-section. However it may be noted that many of the 
blood vessels change their cross-section slowly along their 
lengths [Schneck et al. (197c:)] and may have multiple stenoses 
at junctions, bends etc. Talukder et al. (1976) have cc»iducted 
experiments in the case of blood flow in artery with multiple 
stenoses and pointed out that the pressure drop increases 
linearly with numbers of stenoses. Keeping this in view# in 
this chapter a model for the flow of a power law fluid through 
tubes of non-uniform cross-section with two stenoses is investi- 
gated. The non-uniform distribution of the red cells across 
the tube and the existence of the peripheral layer have beei 
accounted for by considering radial variation of the consistaicy. 
Solutions for mild stenoses have been obtained and the effects 
of heights of stenoses# flow behaviour index etc. on the 
resistance to flow and wall shear stress have been studied. It 
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has been shown that the model presented here, predicts the 
results in conformity with the experimental work of 
Talukder et al* (1976), 

6.2 ANALYSIS : 

Consider the flow of blood (non— Newtonian power law fluid) 
through an artery (circular tube) with two stenoses, the 
non-stenosed parts of which consist of regions with uniform 
radius and slowly varying radius R (z) as slx)wn in Rig. 6.1. 
It is assumed that the stenoses have developed synmetrically and 
are very mild, the shapes of which can be represented by cosine 
curves [ Young (1968) Thus, the varying radius R(z) of the 
tube can be written as follows * 


R(z) 


j^o “ L^ [l+cos£^ (z-dj_-2 )] 

I 

I ° 

) 6 

i ^o - / [H-cosJ2(2_b^)] 

6 2 

R*(z) - ^ [1+cos ^(z-B^)] 


/ O < z < 

; d^+L^ < z < ^ 

L 

' “l 2 ' - ^ ~ 1 

; B^ < z 1 
L 

; b 3_ + ^ i z < B 

( 6 , 1 ) 


Here, L^, and 6^ (i = 1,2) are the lengths and the maximum 
thickness of the two stenoses (the suffixes 1 and 2 here indicate 
the first and second stenosis respectively) and are such that 
the restrictions for mild stenosis [ Young (1968) J are satisfied. 




J6S 


r6j « min (R -R 
f ^ o' out 


L«i 

where R^^^ = R(z) at z = B, 


(i = 1,2) 


( 6 -. 2 ) 


Following Young (1963) and Shukla et al. (1980b) we ccffisider 
the flow to be steady, laminar and axially synmetric. As in 
this analysis blood has beoi characterized by the power— law 
model whose consistency may vary along the radial direction, 
the basic equation governing the flow is 


r I? = 


(6.3) 


where the shear stress t is given by 


= -m(r)l M 1“-^ -I? :=m(r)(- 


< o, R > r > o 


(6.4) 


Here w is the axial velocity, p the pressure, n is power-law 
exponent which is assumed to be constant and may be taken as 
the mean of the flow behaviour index over the radius, m(r) 
is the consistency function which varies radially. 

Combining eqns. (6-. 3) and ( 6 . 4 ) we get our basic equation 


i|p[ria(r)C- 9^)"] + „ 


(6.5) 


The boundary conditions for w are provided by the no slip 
condition at the boundary of the tube and the axial s^^nmetry 
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of the flow 


5w 

dr 


= o 


w = o 


at r = o 
at r.. = R (z) 


( 6 . 6 ) 


Solving equation ( 6 . 5 ) and using the conditions (6.6) the 


axial velocity can be obtained as" 

The volumetric flow rate is defined asF 

R(z) 

Q = / 2 7C rw dr 

o 


(6.7) 


( 6 . 8 ) 


which on substituting w from eqn. (6.7) into eqn. (6,8) is 
given by 


where 


If =-2 

R(z) 2+1/n 

I(z) = / ' y dr 


(6.9) 

( 6 . 10 ) 


Integrating eqn. (6.9) with respect to z we get the pressure 
difference Ap along the total length of the tube as follovjs 

( 6 . 11 ) 


B 


Ap = 2(“)’^ / 


_d z , _ 
o (Kz)]’^ 

The resistance to the flow X is calculated from ( 6 .II) as 


follows 


X _ ^ = 2 n ® . dz 

° Q * o 


( 6 . 12 ) 


The shearing stress x « on the wall of the tube (r = R(z)) 


w 
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is given by [Man ton (l97l)] 

l+CdR/dz)^ ^ r^iz) 

which on using eqns. (6»7) and (6»9) can be obtained as' 
_ l-(dlVdz)^ Td r w Q -.^1 


(6-13) 


(6.14) 


Let and be the resistance to the flow and the shearing 
stress at the wall respectively for the flow of a power law flui<J 
(with fluid behaviotir index n and constant consistency with 
flux in a' tube' of uniform radius and length B having no 
stenosis. These can be obtained from (6.12) and (6.14) 


respectively as 


Q " n ^ ^3n+l 

o 


(6.15) 


(3n+l)n ._1 


(6.16) 


Thus from equations (6.12), (6.14) and (6 .15 )/ (6.16) for same 
flux Q/ the non -dimens i on ali zed resistance and the shear stii^ss 


can be obtained as follows- • 

R3n+1 g 

Sn+l’ B - ’ 




(6.17) 


- ^w 1- (dR/dz) ^ / n \ n( ) fx 

6 .3 Effect of the Peripheral Layer I 


(6.18) 


To study the effect of the peripheral layer consistency on 



K and ir we consider the flow to be consisting of a core region 
with consistency surrounded by a peripheral layer of ccxisist— 
ency m 2 * Thus the consistency function m(r), in this case can 
be written as 

m, 0 r 

m(r) = { (6.19) 

= Pm-, R (z) 1. r <. R(z) 

^ 1 p “ — 

where R,^(z) is the radius of central core and peripheral layer 

Jr 

has thiclcness (R(z) — R (z)). 

Jbr 

Following Shukla et al. (l980c) the central layer radius 
R (z) is given by 

R (z) = cxR(z) (6.20) 

P 

where oc 1 is a constant. 

Using (6.19) and (6*20) and integrating (6.10) we get 

[I(z)]’^ = R^^'*'^(z) (6.21) 

where 

G = (|) [1- (g^2Z) 

Substituting the expression for l(z) from (6.21) into 
(6,17) and (6.18) and introducing the following non— dimension ali- 
zed variables 

( Z/ L^/ B^) = (z/dj^, ^1^/ 

R(z) = R(z)/ Rq,R*(z) = R*(z)/ Rq# 6^ = 6^/Bq Ci =1,2) 



we finally get 


1 

; 

o 


di 


(6.23) 



(6.24) 


l-(dR/dz)^{R /B)^ 

X = - - o . _ , 1.^ 

l+(dj^/di)^(R^/B)^ G[R(2)]^’^ 

6.4 Results and Discussions * 

It can be seen that the results presented in this paper 
reduce to those of Shukla et al. (1980) when the variable r^ion 
is abstent and P = 1. To explicitly see the effects of various 
parameters, the resistance to the flow and the shear stress 
acting on the wall have been numerically computed by taking 

R* Cz)/ R^ = exp[_ K( 2 -B^)^] = exp[lc(2-5^) 
where K = KB^ 

and 

d^ = 0.2, = 0.2, ^2 = ^*2, = 0.8, R^/B = 0.1,a = 0.95 

The effects of heights of stenoses 6^^ and &2 * .power law 
index n, wall exponent paraueter K and consistency ratio P on 
resistance to the flow X are shown in Pigs. (6.2)- (6.4). It- 
can be noted from these figures that the resistance increases 
with increase in heights of stenoses 6^^ and also with 

the increase in consistency ratio P. The resistance ratio 
decreases with decrease in power law index n. It can also be 
observed that for given values of P and n resistance increases 
for a convergent tube (K<0) ccmpared to its values for a 
straight tube (K = O) but it decreases for a divergent tube 
(K >0). It can also be noted that X increases with number 
of stenoses. This observation is in conformity with the 
experimental findings of Talukder et al (1976) . 
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Fig.6'4 Effect of $2 on ^ with variation in 0 and K 
(n = 0.75, (S) = 0.1) 



J 77 



Fig.e-SEff^ct of 62 on T 2 ( n = 0 . 75 , 0 = 0.5 ) 
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Fig .6*6 Ef Icct of B and n on T 2 ( S 2 -0.1 i k = 0 ) 
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The shear stress acting over the length of the second 
stenosis (in the region I < L 2 / 2 ) is shown graf^dcally 

xn Pigs. (6.5)-(6.6) for different values of 62 / n and P. It 
can be observed that the wall shear ratio increases with 
increase in power law index n. The wall shear also increases 
with increase in P. It is seen that the value of f 2 oscillates 
along the length of the stenosis having a maxiitnm value at the 
mid point of stenosis and two mintoa at l/ 4 th and 3 / 4 th of 
the length of the stenosis* Ihe increase in maximum and decrease 
in minimum values of T 2 is very significant as 62 increases* 



Effect 

of K on T 2 ^ 

62 = 0 . 10 , 

P = 0.5') 






n 

ic 


- ' - . 

- , - 



O-.l 

0.3 

0.5 


- 0.1 

0.6542 

0.5604 

0.6524 

0.25 

0.0 

0,654^1 

0.5599 

0.6523 


0.1 

0.6540 

0.5595 

0.6521 


- 0.1 

0.9482 

0.7427 

0.8205 

ID 

* 

0 

0.0 

0.9480 

0.7420 

0.8200 


0.1 

0.9479 

0.7413 

0.8195 


- 0.1 

1,0473 

0.7845 

0.8442 

0 

0 

0,0 

1,0471 

0.7837 

0.8435 


0.1 

1,0469 

0.7829 

0.8429 


The effect of K on T 2 is shown in Table 6,1, for fixed 
values of P,n and 62 * Since the cross sectional radius of 
the tube does not vary upto half of the length of the second 
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stenosis, the parameter K does not affect in that regicm 
and thus its values in that region are not shown in table 6 *1 
It can be noticed that increases for a convergent tube and 
decreases for a divergent tube in comparision to its value 
for a straight tube. 

6,5 CONCLUSION I 

A matonatical model for the steady flow of a ix>n-Newtoniain 
power law fluid with varying cross~section and having tw 
stoaoses, one located at the Junction, has been presented. 
Solutions has been obtained for mild stenoses case and it has 
been shown that the resistance to the flow increases with 
number of stenoses, for a given flux. This result is in 
conformity with the experiment by Talukder et al» (1976), 
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